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ABSTRACT 

We study the relation between the surface density of gas and star formation rate in twenty 
moderately-inclined, bulgeless disk galaxies (Sd-Sdm Hubble types) using CO(l-O) data from the 
IRAM 30 m telescope, H I emission line data from the VLA/EVLA, Ho data from the MDM Obser- 
vatory, and PAH emission data derived from Spitzer IRAC observations. We specifically investigate 
the efficiency of star formation as a function of circular velocity (i> c irc)- Previous work found that the 
vertical dust structure and disk stability of edge-on, bulgeless disk galaxies transition from diffuse dust 
lanes with large scale heights and gravitationally-stable disks at i? c j rc < 120 km s -1 (M* < 1O 1O M0) 
to narrow dust lanes with small scale heights and gravitationally- unstable disks at u c i rc > 120 km s -1 . 
We find no transition in star formation efficiency (Esfr/^hi+h 2 ) & t Vtirc = 120 kms -1 , or at any 
other circular velocity probed by our sample (« c j rc = 46 — 190 kms -1 ). Contrary to previous work, we 
find no transition in disk stability at any circular velocity in our sample. Assuming our sample has 
the same dust structure transition as the edge-on sample, our results demonstrate that scale height 
differences in the cold interstellar medium of bulgeless disk galaxies do not significantly affect the 
molecular fraction or star formation efficiency. This may indicate that star formation is primarily af- 
fected by physical processes that act on smaller scales than the dust scale height, which lends support 
to local star formation models. 

Subject headings: galaxies: spiral — galaxies: star formation — galaxies: ISM — radio lines: galaxies 



1. INTRODUCTION 

High resolution studies of the Milky Way and nearby 
galaxies show that star formation mainly occurs in giant 
molecular clouds (GMCs), with typical GMC masses be- 
tween about 10 3 and 10 7 M and ra dii between about ten 
and several hundred parsecs (e.g., iFukui fc Kawamural 
120101) . It is challenging to study molecular gas on 
these scales in even the nearest galaxies (< 4Mpc; e.g., 
iBolatto et al.l l2008) and studies of more distant galaxies 
can only currently investigate the average properties of 
gas and stars on larger scales. An advantage of more 
distant galaxies is that they exhibit a much wider range 
of physical properties, which enables investigation into 
the influence of environment (e.g., mid-plane pressure, 
metallicity, and scale height) and processes that act on 
large scales (e.g., shear and large-scale gravitational in- 
stabilities) on star formation. 

Extragalactic studies of star formation have found that 
the star formation rate (SFR) surface density (Esfr) 
and gas surface density (S gas ) are correlated in the form 
of the Kennicutt- Schmidt law: Esfr oc (jSchmidtl 
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119591 : iKennicuttl |T998) . This star formation law has 
been studied by averaging over scales as small as about 
100 pc and as large as the entire optical disk. Recent 
high-resolution studies have found that the SFR surface 
density is more strongly correlated with the molecular 
gas surface density (£h 2 ) than with the atomic gas sur- 
face density (Shi), with Esfr oc EjL and N be tween 
0.8 and 1.5 (IWong fc Blitd l2002t Kennicutt et~aT]|2007 



Bigicl et al.ll20 08; Blan c et al 



Kennicutt et al. 
2009: Schruba etal 



2010; 



Liu et all 120111 : iRahman et alJl2011l ). This result con 
firms, from an extragalactic perspective, that stars form 
from molecular gas and has led to an expansion in the 
scope of many star formation studies to investigate how 
environment and large-scale processes affect the molecu- 
la r fraction in the in terstellar medium (ISM). 

iLeroy et al.l (|2008h addressed environmental effects on 
the molecular fraction and star formation efficiency 
(SFE) with high-quality, 750 pc-resolution observations 
of Ssfr, Shi, and £h 2 over the optical disk of 23 nearby 
galaxies. They compared these observations to many star 
formation models and thresholds and concluded that no 
model fit the data sufficiently well to be declared a clear 
favorite. This result led them to suggest that physics 
that acts on scales smaller than their resolution is most 
important for setting the molecular fraction and SFE. 

While no model w a s an ideal fit to the data pre- 
sented in iLerov et al.l ()2008| ). the best fit was arguably 
with a model in which the ratio of molecular to atomic 
surface density is related to the mid-plane pressure 
(Ph): R mo i = Eh 2 /Shi oc P£, and the molecular SFE 
(SFE[H 2 ] ee S sfr /Eh 2 ) is constant , such that S FE = 
Esfr/Ehi+h 2 = SFE(H 2 ) ^"ji - ILerov et al.1 (|2008h 



noted that a model where the molecular SFE is constant 
requires the population of GMCs in any region to be 
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sampled from the same distribution of properties (e.g., 
size and mass), independent of environment. Further- 
more, once GMCs form, the general environment cannot 
have a strong influence on their properties. Finally, one 
must compare the model to observations with a number 
of GMCs per resoluti on element so as to average over 
evolutionary effects. lEhnegreenl (|1993l ) predicted that 
i?moi should depend on the mid-plane pressure and the 
interstellar radiation field (J): i? mo i oc P^' 2 j -1 ■ As- 
suming Ssfr °c ^h 2 and j oc Ssfr. the mod e l pre - 
dicts P mol oc P h Q , with a = 1.2. IWong fc Blit3 ()2002f ) 
studied seven molecu lar-dominated spiral galax ies and 
found a = 0.8, while iBlitz fc Rosolowskvi |2006f ) found 
a = 0.92 in their study of fourteen galaxies with a large 
range of i? mo i- V alues between a = 0.5 and 1.2 encom- 
pass most of the iLerov et"al"1 (|2008l) data on the SFE 
v ersus Ph plane. 

lOstriker et al.l (|2010f) recently presented a star forma- 
tion model that produces an approximately linear re- 
lationship between P mo i and Ph- The authors divided 
the ISM into a diffuse component and a gravitationally 
bound component. The fraction of gas in each compo- 
nent is set by the requirement that gas pressure in the 
diffuse component is balanced by the gravity of stars, 
dark matter, and gas (both diffuse and bound), while 
heating (mainly by ultraviolet (UV) photons from O and 
B stars formed in the bound component) balances cool- 
ing. The model assumes that the SFE wi thin the grav- 
itation ally bound component is constant. iBolatto et al.l 
(120 111) added a metallicity-dependent heating term to the 
lOstriker et al.l (|2010l) model, which resulted in agreement 
between the large H I surface densities observed in the 
Small Magellanic Cloud and the su rface density of diffuse 
gas calculated with the adjusted lOstriker et al.l (|2010f ) 
model. 

Another lead i ng st ar formation model is that of 
iKrumholz et al.l ((2009), where the molecular fraction is 
determined by processes that act on scales no larger than 
~ 100 pc, which is the size of atomic-molecular complexes 
in the model. Specifically, the molecular fraction is set by 
the balance between the formation of molecular hydro- 
gen on the surfaces of dust grains and the destruction of 
molecular hydrogen by UV photons. Both dust shielding 
and H2 self-shielding contribute to the survival of molec- 
ular hydrogen in the interior of the complexes. Stars 
form only from molecular gas and the model produces a 
constant molecular SFE because properties of molecular 
clouds, like Eh 2 , are independent of general ISM condi- 
tions, at least while gas surface densities in the general 
ISM are less than CMC densities (~ 85M Q pc- 2 ). 

In this paper, we study star formation in a sample 
of 20 bulgeless disk galaxies. Bulgeless galaxies are in- 
teresting from a number of perspectives. Their exis- 
ten ce in relatively large numbers (>15% of d isk galax- 
ies; [Kautich^rialJUool; |Kormendxii!]d][2Q]^ provides 
an important constraint on hierarchical galaxy formation 
models, in which galaxies generally have rich merger his- 
tories that lead to bulge growth. The agreement between 
models and observations is becoming better as feedback, 
high gas fractions, and satellite mergers on radial or- 
bits are included in the models (Robertson et al. 2006, 
Hopkins et al. 2008, Brook et al. 2011; but see also 
Scannapieco et al. 2011). Under the assumption that 
bulges do form when significant merger events occur, 



bulgeless galaxies are a suitable sample to study secular 
evolution, where internal processes like star formation 
lead to changes in the galaxi es, such as bulge growth 
(|Kormendv fe Kennicuttll2004l ) 

A number of works have studied the co mponents of 
star f ormation in late- type disk galaxies. iBoker et al.1 
(2003) found that their sample of 47 late-type spirals 
are similar to earlier-type spirals in that they fall on 
the approximately linear correlation between far-infrared 
(FIR) luminosity, which traces star formation, and the 
molecular hy drogen mass wit h in the central few kpc. 
Furthermore, iMatthews et al.1 (|2005h found that low- 
surface brightness, late - type disks lie on this same re- 
lation. iDalcanton et al.l (|2004f ) studied the dust and cold 
ISM structure in a sample of 49 edge-on bulgeless disk 
galaxies. They inferred that there is a sharp transition 
in dust lane structure with circular velocity (w C irc) based 
on measurements of R — K color versus height above 
the midplane. Galaxies with u C j rc < 120 km s -1 (we also 
refer to these as low-i> c i rc galaxies) appear to have no 
dust lanes while galaxies with v c i IC > 120kms _1 (high- 
^circ galaxies) have well-defined dust lanes (the stellar 
mass Tully-Fish er relation relates v c \rc = 120 km s" 1 with 
M* - 10 10 M g ; lBell fc de Jong||200lD . The authors con- 
cluded that the transition is likely due to a transition in 
dust scale height rather than due to a sharp transition 
in the quantity of dust present: low-w C j rc galaxies have 
diffuse dust lanes with large scale heights while high-t; c i rc 
galaxies have dust lanes with smaller scale heights. They 
came to this conclusion because the dust structure tran- 
sition occurs over a relatively narrow range in circular 
velocity (~ lOkms -1 ), and therefore over a relatively 
narrow range in gas and total mass, where the dust-to- 
g as ratio (PGR) d o es no t vary substantially. 

IDalcanton et al.l (|2004l ) also found that disk stability, 
parametrized by a generaliz ed Toomre Q parameter in- 
cluding both gas and stars (|Rafikovll2001| ). is correlated 
with the dust structure, with low-i> c i rc galaxies generally 
stable and high-?; c i rc galaxies generally unstable. Fur- 
thermore, they concluded that a sharp change in the con- 
tribution of turbulence to the stability parameter is the 
likely cause of the stability and dust scale height tran- 
sitions. The authors suggested that high-w c i rc galaxies 
have turbulence dominated by supernovae explosions and 
gravitational instabilities while low-w C i rc galaxies have 
turbulence dominated by only supernovae. Independent 
of the source of the turbulence, the turbulent velocities 
must be lower in the high-u c j rc galaxies to explain the 
stability results. An alternative interpretation for the 
dust structure transition is that it is due to differences 
in stellar surface density and dust opacity, as suggested 
b vlHunter fc Elmegreenl (120061) . 



Dalcanton et al. (20041) noted that the cold, star- 



forming gas should have a similar distribution to the 
dust, with larger scale heights in low--u C j rc galaxies com- 
pared to high-w c j rc galaxies. They hypothesized that a 
transition in SFE might accompany the transition in dust 
scale height and stability if the volume density of gas is 
the relevant quantity for setting the SFR. A low-w c i rc 
galaxy with a larger scale height but the same £ gas rel- 
ative to a high-i> c i rc galaxy will have a lower gas volume 
density (p gas ). The low-i; c i rc galaxy likely also has a lower 
gas pressure because pressure is proportional to the gas 
volume density (P oc p gas Cg as , where <7 gas is the gas ve- 
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locity dispersion). In the context of the star formation 
model where the molecular fraction is set by the mid- 
plane pressure, we then expect lower R mo \, Ssfr, and 
SFE in the low-w c ; rc galaxy. 

In this paper, we address whether there is a SFE tran- 
sition at w c ; rc = 120 km s -1 in our sample of bulgeless 
disk galaxies ( Section 14.2)) . We also investigate whether 
scale height differences affect the SFE and discuss impli- 
cations for the scale of physical processes that are im- 
portant for setting the molecular fraction and the SFE. 
We examine our results in light of recent star formation 
models such as the the mid-p l ane p ressure model and 
the model of iKrumholz et all (|2009f ) (Sections |4~3] and 
14. 4p . To carry out this study, we trace molecular gas with 
CO(l-O) data from the Institut de Radioastronomie Mil- 
limetrique (IRAM) 30 m telescope and atomic gas with 
H I 21 cm data from the Very Large Array (VLA; Wat- 
son et al. 2011, hereafter Paper I). We trace the SFR 
with Ha data from the 2.4 m Hiltner Telescope of the 
MDM Observatory combined with polycyclic aromatic 
hydrocarbon (PAH) emission data derived from Spitzer 
Space Telescope Infrared Array Camera (IRAC) observa- 
tions. We also estimate the stellar mass from the Spitzer 
IRAC data. These observations, and measurements de- 
rived from the data, are described in Section [2j We de- 
scribe the quantities that we derive from these measure- 
ments in Section [3] Our results are in Section |4] and we 
discuss these results in Section [5] 

2. SAMPLE SUMMARY, OBSERVATIONS, DATA 
REDUCTIONS, AND MEASUREMENTS 

Our sample is composed of 20 Sd-Sdm galaxies 
within 32 Mpc, with circular velocities between 46 and 
190 km s" 1 . Th e se pro perties are well matched to the 
IDalcanton et al.l (120041) sam ple. However, in contrast to 
the iDalcanton et al.l ( 20041 ) sample of edge-on galaxies, 
we selected our galaxies to be moderately inclined, with 
inclinations between 16° and 56°, such that we can ac- 
curately measure the SFR and gas surface densities and 
place the galaxies on the star formation law. Section 2 in 
Paper I and Table [T] provide a description of the sample 
selection. 

Many of the measurements described in this section 
were carried out to derive surface densities - of gas, SFR, 
and stars. These surface densities must be measured over 
the same area. The IRAM CO(1-0) data are the limiting 
factor, as they are single-beam, FWHM = 21" measure- 
ments centered on each galaxy. Therefore, we measured 
the emission within a 21"-diameter circular aperture cen- 
tered on the IRAM pointing center, the coordinates of 
which are listed in Table [2j for the following datasets: 
the H I data from the VLA, the Ha data from the MDM 
Observatory, and the PAH and 4.5 /im data from Spitzer 
IRAC. 

2.1. IRAM 30 m CO(1-0) 

Thirteen of our objects were observed in the CO(1-0) 
and CO(2-l) lines at 115 and 230 GHz in May 2007 with 
the IRAM 30 m telescope on Pico Veleta. Dual polariza- 
tion receivers were used at both frequencies with the 512 
x 1 MHz filterbanks on the CO(1-0) line and the 256 
x 4 MHz filterbanks on the CO(2-l) line. The observa- 
tions were carried out in wobbler switching mode with 
a wobbler throw of 200" in the azimuthal direction. At 



the beginning of the observations the frequency tuning 
was checked by observing a bright galaxy at a similar 
redshift. Observations of the same calibration source on 
different days allowed us to check the relative calibra- 
tion, which was excellent (better than 10%) for CO(1-0). 
The calibration in CO(2-l) was equally good, except for 
one day when the calibration observation was different 
by ~ 35%. Pointing was monitored on nearby quasars, 
Mars, or Jupiter every 60 - 90 minutes. During the ob- 
servation period, the weather conditions were generally 
good, with pointing better than 3". The typical system 
temperature was 300-500 K at 115 GHz and 500-1000K 
at 230GHz on the scale. At 115 GHz (230 GHz), 
the IRAM forward efficiency, F cS , was 0.95 (0.91), the 
beam efficiency, B c g , was 0.75 (0.54), and the half-power 
beam size is 21" (11"). All CO spectra and line intensi- 
ties are presented on the main beam temperature scale 
(T m b) which is defined as T mb = (F e g/B e g) x T£. For the 
data reduction, we selected the observations with good 
quality (taken during satisfactory weather conditions and 
showing a flat baseline), averaged the spectra from the 
individual scans of the source, and subtracted a constant 
continuum for the CO(1-0) spectra and a linear contin- 
uum for the CO(2-l) spectra. 

Figure [T] shows the CO spectra for the objects observed 
in May 2007. The black spectrum shows the CO(1-0) 
data and the red spectrum shows the CO(2-l) data. We 
did not use the CO(2-l) data in this paper, but show it 
for completeness and to corroborate some of the weaker 
CO(1-0) detections. The solid black horizontal line is 
the velocity range over which we integrated to derive the 
CO(1-0) line intensity. The dashed horizontal line is 
centered on the systemic velocity of the galaxy, derived 
from velocity field modeling of our H I emission line data. 
The width of the dashed line is the width of the H I 
emission line at 20% of the peak flux density (W20; see 
Table [1] for these values). The integrated CO(1-0) line 
intensities are presented in Tabled 

We did not detect ESO 544-G03, ESO 418-G008, 
ESO 501-G023, or UGC 6446. For these galaxies, we 
quoted the CO line intensities as upper limits, computed 
with I co < 3 CT rms (W20 Aw) 1 / 2 , where Crms is the noise in 
the spectrum in K and Av is the CO (1-0) spectrum chan- 
nel width, which is lOkrns -1 in the Hanning-smoothed 
spectra of Figure [TJ 

CO(1-0) and CO(2-l) spectra and line intensities for 
six of our objects - PGC 3853, NGC 2805, NGC 3906, 
N GC 4519, NGC 596 4, and NGC 6509 - were presented 
in iBoker et alJ (|2003t ) , also based on IRAM 30 m obser- 
vations and reduced in the same manner. We did not 
obtain CO data for ESO 555-G027. 

2.2. VLA HI 
2.2.1. Integrated H I Line Intensity 

The H I 21 cm data for our sample were obtained from 
the VLA/Expanded VLA, operated by the National Ra- 
dio Astronomy Observatory 8 , for projects AZ0133 (car- 
ried out in 2001 August), AL0575 (carried out in 2002 
June and November), AM0873 (carried out in 2006 Oc- 
tober and November), and AM0942 (carried out in 2008 

8 The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 
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Fig. 1.— CO(l-O) (black) and CO(2-l) (red) spectra for the objects observed in May 2007. See IBoker et al.l 112003 1 for the remaining 
spectra. The black horizontal line designates the velocity range over which we integrated to derive the CO(l— 0) line intensity. ESO 544-G03, 
ESO 418-G008, ESO 501-G023, and UGC 6446 do not have this line because they were undetected. The dashed horizontal line is centered 
on the systemic velocity derived from velocity field modeling of our H I data, and has a width of W20 > which we derived from the integrated 
H I line profile. 



May and 2009 July and August). The galaxies were ob- 
served in the C or CnB configurations, which provide a 
nominal angular resolution of 13". The channel width 
is generally 5.2 km s . The observations and reductions 
are described further in Paper I. We measured the in- 
tegrated H I line intensity from velocity-integrated in- 
tensity maps (with units of Jy beam -1 kms -1 ) created 
from naturally-weighted data cubes. The beam major 
axis FWHM (B maj ) and minor axis FWHM (B mia ) for 
each data cube are listed in Table 3 of Paper I. 

Ten of our objects have i? ma j and B m i n < 21". To 
adjust the H I beam to match the CO(1-0) beam, we 



used the AIPS task CONVL to convolve the integrated 
intensity map so the output Gaussian beam has B ma] = 
B m in = 21". To approximately match the H I beam to 
the CO (1-0) beam for the remaining ten objects with 
£ ma j > 21" (on average, 5 maj = 29") and B min < 21", 
we either used the original integrated intensity map or 
convolved the map to have B m \ n — 21" and the beam 
major axis approximately equal to the original -B ma j. 

Using the resulting maps, we measured the average 
H I line intensity in Jy beam -1 kms -1 within the 21"- 
diameter circular aperture ((/hi)) with the AIPS task 
IMEAN. These values are listed in Table [2j The inte- 
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grated H I line intensity in K km s 1 is given by: 



(Jhi) [Kkms- 1 ] 



6.07 x 10 5 (/hi) [Jybeam^kms- 1 ] 



-^maj ^min 



where B ma j and B 



(i) 

mm are in arcseconds and now refer to 
the beam of the map on which we made the measure- 
ments. These values are listed in Table [2j 

For the objects with B ma j > 21", we assumed that 
the average line intensity within a 21" beam is equal to 
the measured line intensity from our image with a larger 
beam. We estimated the uncertainty introduced by this 
assumption by convolving the integrated intensity maps 
of three objects with B ma j < 21" such that the convolved 
beams match those of the ten objects with B ma j > 21". 
We measured the average line intensity within the 21"- 
diameter circular aperture and compared this to the true 
value measured from the map with a 21" beam. We 
found that the line intensities in our test cases differ from 
the true values by up to 11% (in the case of the NGC 4519 
beam, which has a B ma j = 51.91") and by 4% on average. 
The test measurements both over and under estimate 
the true value depending on the emission distribution. 
Therefore, we included this in our uncertainty estimate, 
but make no correction. 

The main contributors to the final uncertainty in our 
H I line intensities are flux calibration (5%), aliasing (up 
to 11% and described in Paper I), and using an image 
with a beam larger than 21" to estimate the line intensity 
within 21" (on average 4%). Not all objects are subject 
to the latter two uncertainties. Nonetheless, we conserva- 
tively assigned the quadrature sum of these uncertainties 
(13%) as the generic uncertainty associated with our H I 
line intensities within the 21"-diamcter aperture. 

NGC 6509 shows H I in absorption on the east side of 
the galaxy because it is in the foreground of the radio 
source 4C +06.63. The average line intensity within 21" 
is unaffected because the eastern edge of the aperture 
and the western edge of the radio lobe are separated by 
about 20". 



2.2.2. Epicyclic Frequency 

We calculated a representative epicyclic frequency (k) 
for the 21"-diameter circular aperture to use in the stabil- 
ity analysis of Section l4~3l We used K = \/2(l + (3) (v/r), 
where j3 — rflog(u)/rflog(r). We determined j3 and 
the rotation velocity, v, at a radius, r, of 5.25" (half 
the radius of the 21"-diameter aperture) from the H I 
rotation curves presented in Paper I, where we fit a 
tilted ring model to the data to derive v at radii every 
(S maj B min ) - 5 , beginning at (5 majJ B min ) a5 /2. We sim- 
ply fit a line between the origin and the first point in the 
rotation curve, the average radius of which is 8.3". We 
used the slope of the line as an estimate of dv/dr and 
evaluated v at 5.25". There may be inaccuracies intro- 
duced to k evaluated in this manner because the origin 
and first point in the rotation curve are within a single 
beam. Therefore, we also calculated the epicyclic fre- 
quency at 10.5" using the same method as above, but by 
fitting the line between the first and second points from 
the rotation curve, where the average second radius is 
25.0". The epicyclic frequencies calculated at 10.5" are 
smaller than the values at 5.25" by 30% on average. We 
used the epicyclic frequencies evaluated at 5.25" in our 
stability analysis (and list these in Table H]), but include 
30% uncertainties on the values. Beam smearing, where 
many velocity components are within the spatial beam, 
is likely in effect in this region. Beam sm earing leads to 
underestimated velocities and gradients ( S w aters! 1 1999f ) 
and therefore k may also be underestimated. To account 
for this, we use n evaluated at 5.25", as these values are 
larger. 

2.3. MDM Ha 

Ha and continuum images of the galaxies were ob- 
tained at the MDM 2.4m Hiltner telescope over the 
course of four observing runs in January 2007, November 
2007, May/June 2007, and January 2008. Each galaxy 
was observed for between 30 min and 2.5 hours through 
a pair of matched, custom-made 15 nm wide narrow- 
band filters centered at 663 nm and 693 nm, hereafter the 
663bpl5 and 693bpl5 filters, respectively. The Ha emis- 
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sion line falls within the 663bpl5 bandpass for all of the 
galaxies in this sample. Observations were obtained with 
the Direct CCD camera "Echelle," which has 2048x2048 
pixels. The CCD was binned over 2x2 pixels to produce 
a plate scale of 0.55" /pixel, which was well-matched to 
the 1" — 1.5" image quality measured on most nights. 
The field of view was 9.4'. Conditions were photometric 
for most of the observations and a series of spectropho- 
tometric standards were observed for flux calibration, as 
were a series of twilight flats. The exposure times and 
observation dates are listed in Table [3l 

Overscan subtraction, flat fielding, cosmic ray rejec- 
tion, and bad pixel removal were performed with IRAF 9 . 
All of the 663bpl5 and 693bpl5 images of each galaxy 
were regi stered to a common coordinate system with 
SCAMP (jBertinl [20061 and a combined image for eac h 
filter was constructed with SWARP (jBertin et al.ll2002ft . 

Observations of spectrophotometric standard stars 
were used to determine the absolute flux calibration for 
the 663bpl5 filter (including an airmass correction) and 
the relative throughput of the two filters. We calculated 
the expected ratio of stellar flux in these filters for late- 
type galaxies by using the SYNPHOT package in IRAF 
to convolve the filter transmission f unctions with a se- 
ries o f Bruzual & Chariot 1995 (e.g.. iBruzual fc Charlotl 
2003) stellar population synthesis models with both con- 
tinuous and exponentially declining (t = 1 Gyr) star 
formation histories. The range of flux ratios from these 
models and the relative throughput of the two filters were 
used to scale the 693bpl5 image and subtract it from the 
663bpl5 image to create an Ha image for each galaxy. 

We measured the Ha flux for each galaxy within the 
21"-diameter circular aperture and within a circle of di- 
ameter -D257 t ne B-band major isophotal diameter at 
25 magarcsec -2 , with aperture photometry. All Ha flux 
measurements were multiplied by a factor of 0.75 to 
account for emissi on from [N II] AA6548, 6584 in the 
663npl5 bandpass (|Kennicuttl ll983). The Ha flux mea- 
surements were also correct ed for Galactic ex tinction 
using the extinction law of lO'Donnelll (11994ft assum- 
ing Ry =3.1 and reddening values from Schlcgel et al. 
(1998) and tabulated on NED 10 . The final values within 
the 21" aperture are listed in Table [2] 

A number of our galaxies show substantial variation in 
the Ha flux within the 21" aperture if we vary the center 
of the aperture by the ~ 2 — 3" pointing accuracy of the 
IRAM 30 m. Including this uncertainty, as well as un- 
certainty in the [Nil] correction and absolute calibration 
uncertainties, we estimate the uncertainty in the Ha flux 
within the 21"-diameter aperture to be about 30%. 

Figure [2] compares our Ha fluxes measured within 
D25 r ela tive to publis h ed va lue s from iKoopmann et al. 
(2001), I James et all (12004ft ■ iMoustakas fc Kennicutt 
(2006(), and lEpinat et alJ (|2008). where we corrected the 
published values for [N II] emission and Galactic extinc- 
tion if necessary. Our values are on average 20% smaller 

9 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 

10 The NASA/IPAC Extragalactic Database (NED) is operated 
by the Jet Propulsion Laboratory, California Institute of Tech- 
nology, under contract with the National Aeronautics and Space 
Administration. 
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Fig. 2. — Comparison between our Ha fluxes measured within 
a circle of diameter D25 and published Ha fluxes. Both values are 
corrected for [N II] emission and Galactic extinction. The dotted 
line shows equality. The largest outlier is NGC 0337. 

than the published values. This systematic offset is com- 
parable to our calibration uncertainty and may be due to 
differences in the filter profiles, the angular extent of the 
apertures, and the treatment of bright stars with poten- 
tially large subtraction residuals. We note that this offset 
corresponds to an uncertainty in the star formation rate 
surface density of less than 0.1 dex, which is smaller than 
the 0.2 dex uncertainty that we assign due to the variable 
contribution to dust heating from non-star forming pop- 
ulations. 

2.4. Spitzer IRAC 

Fourteen of the 20 galaxies were observed as part 
of our Spitzer Cycle 5 Program 50102; the remaining 
six (NGC 337, UGC 1862, ESO 418-G008, NGC 2805, 
NGC 4519, and NGC 4561) were observed for various 
other programs. Observations for Program 50102 con- 
sisted of five, dithered observations with a frame time of 
30s in all four IRAC channels (3.6 /im, 4.5 fim, 5.8 /im, 
and 8 /im) . Data for the other six archival datasets were 
comparable to our observations, with the exception that 
the data for UGC 1862 only include the galaxy in chan- 
nels 2 and 4. In the analysis described below, we use the 
channel 1, 2 and 4 data. 

The Basic Calibrated Data (BCD) for all 20 galax- 
ies were processed with Sean Carey's artifact mitigation 
software 11 , which corrects for a variety of effects such 
as muxbleed, column pulldown/pullup, electronic band- 
ing, and first frame effect. We then created mosaics for 
each channel with the Spitzer Science Center's MOPEX 
(MOasicker and Point source Extractor) package. This 
package corrects individual images for background vari- 
ations and optical distortions, and then projects them 
onto an output mosaic image for each channel. These 

11 http:/ /spider. ipac.caltech.edu/staff/carey/irac_artifacts/ 
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mosaic images were used for five measurements for each 
galaxy: the inclination (i), position angle of the major 
axis (PA), the PAH flux density, the 4.5 /jm flux density, 
and the exponential disk scale length. 

2.4.1. Inclination and Position Angle 

We estimated the inclination and PA for each 
galaxy with the IRAF ellipse task, which fits el- 
lip tical isophote s with the iterative method described 
by Uedrzeiewskil (fl987). The 3.6 /jm data (4.5 /im for 
UGC 1862) were fit because this channel has the greatest 
sensitivity to the old stellar population, yet is relatively 
insensitive to dust. 

We found that the isophote fits, particularly in the 
outer, low-surface brightness regions, were relatively sen- 
sitive to the presence of bright stars. We therefore cre- 
ated masks of these star s with the SExtractor package 
(jBertin fc Arnoutslll996l) and included this mask as an 
input to the ellipse task. These masks also excluded re- 
gions with relatively poor coverage in the IRAC mosaic. 
We averaged the PA and ellipticity values for the largest 
isophotes to derive our final PAs and inclinations. These 
values were used as inputs for the rotation curve analysis 
described in Paper I and were reported in Table 5 of that 
work. 

2.4.2. PAH Flux Density 

We calculated the PAH flux density within the 21"- 
diameter circular aperture using the 8 /im images, which 
are dominated by the 7.7 and 8.6 /im PAH features, 
and the stellar emission-dominated 3.6 /im images. We 
measured the 8 /jm and 3.6 /zm flux densities using the 
IRAF task phot. In both measurements, we subtracted 
from each pixel the median sky background, which we 
measured within a large-radius annulus centered on the 
galaxy. We applied the band-specific extended source 
aperture correction to the 8 /zm and 3.6 /zm flux densities 
(described in the IRAC Instrument Handbook; for the 
21"-diameter circular aperture, the correction is 0.985 
and 0.896 for the 3.6 /im and 8 /im bands, respectively). 
Finally, the PAH flux density is the 8 /im flux density, 
less the stellar emission contribution, which we estimated 
by scaling the 3 . 6 nm flux density by 0.255 (as used in 
iKennicutt etlfl1l2009t ). For UGC 1862, we used the 
4.5 /im flux density, scaled by a factor of 0.389, to re- 
move the stellar contribution to the 8 /im band. This 
scale f actor is derived from values quoted in lHelou et al.l 
(2004), except we assumed that all the 4.5 /im-band emis- 
sion is stellar. The uncertainties on the calibrated 3.6 /zm 
and 8 /im flux densities are about 10% (jHora et al.ll200l . 
We quote PAH flux density uncertainties that are sim- 
ple error propagated values. Our PAH flux densities are 
presented in Table [2j 

2.4.3. 4.5 /zm Flux Density 

We calculated the total 4.5 /zm flux density and the flux 
density within the 21"-diameter circular aperture such 
that we can derive the total stellar mass and the stel- 
lar mass surface density in Section |3~51 We first masked 
out bright foreground stars with the IRAF APPHOT 
and DAOPHOT packages, in particular the phot and 
substar tasks. We used large 4.5 /zm to 8.0 /zm flux ra- 
tios to identify foreground stars. For a few bright stars 



that were not adequately subtracted, we manually re- 
placed the affected pixels with values from a region at a 
similar radius from the center of the galaxy. 

We used the IRAF task phot to measure the flux den- 
sity within the 21"-diameter circular aperture. For the 
total flux density, we used the IRAF task ellipse to 
measure the flux density within an ellipse where the semi- 
major axis (SMA) is -D25/2 and the PA and ellipticity 
(e = 1 — cos i) were derived from a combination of rota- 
tion curve analyses on the VLA H I data and the ellipse 
fits of the IRAC 3.6 /zm data, described in Section [2.4. II 
(PA and i are given in Tabled]). In both measurements, 
we subtracted the median sky background from each 
pixel, measured within a large-radius annulus centered 
on the galaxy. 

To confirm that the D25 aperture is appropriate for the 
total flux density measurement, we identified the SMA 
where the surface brightness profile flattens (specifically, 
where the surface brightness decreases by, on average, 
less than 2% over a number of apertures) . The flux den- 
sity enclosed within this ellipse was generally within 7% 
of the flux density within the D25 ellipse, except in two 
cases: UGC 6446 and UGC 6930, where the flux den- 
sity within the D25 ellipse was larger by a factor of 1.9 
and 1.5, respectively. We henceforth use the flux density 
within the -D25 ellipse as the total flux density. 

We applied the extended source aperture correction 
to the 21" and D25 flux density measurements, with a 
value of 1.013 for the 21" aperture and an average value 
of 0.944 for the D25 measurements. We did not apply a 
color correction. The 21" and total 4.5 /zm flux densities 
are listed in Table [2j 

2.4.4. Exponential Disk Scale Length 

We estimated the exponential scale length of each 
galaxy using the IRAF ellipse task on the 4.5 /zm data, 
allowing the center, PA, and ellipticity to vary as a 
function of semi-major axis. We fit an exponential to 
the mean isophotal intensity profile to derive the cen- 
tral surface brightness and scale length. We excluded 
PSF and/or bar components from the profile fit based 
on visual examination of t he images and a provisional 
GALFIT (jPeng et al.ll2002D analysis. We created a two- 
dimensional image representing the ellipse profile with 
bmodel, subtracted the model from the original image, 
and found that the standard deviation in the region of 
the galaxy within the residual image is typically less than 
about ten times the standard deviation in a galaxy-free 
region of the original image. Given the small-scale struc- 
tures present in most of the galaxies, we accepted these 
values and fits. The scale lengths are listed in Table |4] 
and are between 0.69 and 3.4 kpc. 

We assigned an error of 20% to the scale lengths, based 
on comparing the scale lengths computed as described 
above to scale lengths computed from ellipse runs where 
we held the PA and ellipticity fixed as a function of semi- 
major axis at the values from Table[TJ We confirmed that 
scale lengths derived from the 3.6 /zm data are generally 
consistent with the values derived from the 4.5 /zm data 
within the uncertainty (they differ by 9% on average). 

3. DERIVED QUANTITIES 

For all the surface density calculations below, we used 
the 21"-diameter aperture, to match the beam of our 



CO(l-O) data from IRAM. Table H lists the physical size 
of 21" (0.7 - 3.2 kpc) and the parameters derived in the 
following sections. The surface densities are all within 
the deprojected area of the aperture, which we calculated 
with the inclinations from Table [TJ 

3.1. Atomic, Molecular, and Total Hydrogen Surface 

Density 

The H I surface density is given by: 

E m [M Q pc- 2 ] = 0.015 (/hi) cos(i), (2) 

where (/hi) is the average integrated H I line intensity 
within the 21"-diameter circular aperture in Kkms -1 
from Section 12.2.11 We did not include a correction for 
He. The Shi uncertainty is dominated by the contribu- 
tion from the line intensity uncertainty and the typical 
uncertainty in log (Shi) is 0.06 dex. The H 2 surface den- 
sity is given by: 



E H2 [M q pc" 



3.2 



Aco 



2.0 x 10 20 cm- 2 (Kkms 



^— T Ico cos(i), (3) 



where Aco is the CO-to-H2 conversion factor and Ico is 
the CO line intensity in Kkms -1 from Section POI We 

used a constant Aco of 2.8 x 10 20 cm~ 2 (Kkms -1 ) . 
We used this value rather than the Milky Way value 
of 2.0 x 10 20 cm " 2 (Kkms -1 ) 1 so we can plot the 
iKennicuttl ()1998[ ) total hydrogen star formation law rel- 
ative to our data. Again, we did not include a correction 
for He. 

The typical uncertainty in log(£H 2 ) is 0.07 dex, due 
mainly to the CO line intensity un certainty. We did 
not include uncertainty due to Aco- iLerov et"aT1 (|2011l ) 
studied five local group galaxies and concluded that 
Xco is relatively constant at the solar value for 12 + 
log (O/H) > 8.4 and increases with decreasing oxygen 
abundance below 12 + log (O/H) ~ 8.2 - 8.4. We have 
only one galaxy where we estimated that the oxygen 
abundance is below 12 + log (O/H) = 8.4 (Section I3T5]) . 
so we do not expect much Aco variation in our sample. 

We also use the total hydrogen surface density, 
Shi+h 2 = ^hi + Sh 2 • The typical uncertainty in 
log(SHi+H 2 ) is 0.05 dex. 

3.2. Star Formation Rate Surface Density 

We used Ha emission to trace the unobscured SFR and 
PAH emission to trace the obscured SFR. The Ha emis- 
sion is due to recombination in H II regions, which are 
ionized by O and early B stars. PAHs are small dust 
grains (or large mol ecules) that are primarily excited b y 
single UV photons ()Sellgrenlll98llLeger & Pugetlll984D . 
Because ionizing radiation is not required to heat PAHs 
(and in fact th ere is evidence that ionizing radiatio n 
destroys PAHs; iHelou et all 12004 iPeeters et all 120041 ). 
PAH emission traces lower-mass, longer-lived stars than 
those ultimately responsible for the Ha emission. PAH 
emission has b een used to trace the total SFR (e.g., 
IZhu et all [20081) . but there is variation in the 8 /im lu- 
minosity at a giv en SFR due to envir o nment, especially 
metallicity (e.g., ICalzetti et al.ll2007i ). iKennicutt et al.l 
(2009) used 75 galaxies to calibrate SFR estimates based 
on Ha and PAH emission by comparing the combined 



Ha and PAH luminosity to the Ha luminosity corrected 
for extinction using the Balmer decrement. The authors 
found that SFRs calculated with Ha and PAH emission 
agreed with their reference SFRs with as little scatter as 
SFRs calculated with Ha and 24 ^m emission We used 
the SFR calibration of IKennicutt et all (|2009l ) to calcu- 
late the SFR surface density within the 21"-diameter cir- 
cular aperture: 

S S FR[M Q yr -1 kpc -2 ] = 

7.30 x lO^OFka + 1.1 x 10~ 28 v F PAH ) cos(i), (4) 

where Fjj a and /pah are the Ha flux in ergs -1 cm~ 2 
and the PAH flux density in mJy within the 21"-diameter 
circular aperture, v is the central frequency of the 8 /im 
IRAC band in Hz, and the constant includes the aperture 
area. This assumes t he initial mass function (IMF) from 
ICalzetti et al.l (120071) , which is similar to that presented 
in lKroupal(|200lD . The IKennicutt etafl (|2009h SFR was 
calibrated for galaxy-averaged data, but should be ap- 
propriate for our data because we generally probe regions 
that are a couple of square kpc in area and should there- 
fore contain a number of star forming regions. Even using 
both Ha and PAHs to trace the SFR, galaxies with low 
metallicity could have underestimated SFRs because the 
fraction of dust mass in PAHs decreases at low metal- 
licity, particularly below go/4 (12 + log (O/H) < 8.3; 
iDraine et all 120071: ISmith et al1l2007l) . In the following 
section we estimate the oxygen abundance of each galaxy 
from the stellar mass and find that only one galaxy 
(ESO 501-G023) has 12 + log (O/H) < 8.3. Therefore, 
we do not expect significant SFR underestimates in our 
s ample. 

IKennicutt et al.l (j2009h discussed that the dominant 
source of uncertainty in their SFRs is due to varying con- 
tributions to dust heating from older stellar populations 
(> 100 Myr). Based on their suggestions, and including 
the fact that our sample has a limited range in morphol- 
ogy and therefore star formation history, we assigned a 
general uncertainty of 0.2 dex to our SFR surface den- 
sities. This dominates over the contribution due to Ha 
flux and PAH flux density measurement uncertainties. 

In upcoming sections, we use the SFE defined as 
£sfr/£hi+h 2 in yr . The typical uncertainty in the 
SFE is 0.2 dex. 

3.3. Stellar Mass Surface Density, Total Stellar Mass, 
and Oxygen Abundance 

We derived the stellar mass surface density within the 
21"-diameter circular aperture and the total stellar mass 
from the 4.5 /im flux densities. We chose to use the 
4.5 /im data over the 3.6 /im data because there are no 
PAH emission features in the 4.5 /im band. 

To estimate the stellar mass, we used a relation- 
ship between K - band mass-to-light ratio and color from 
iBell fc de Jongj (|2001h . derived from stellar population 
synthesis modeling: 

log (?f) = lA3(J-K s )- 1.53, (5) 

where is the mass-to-light ratio in the K band in 
Mq/Lk,©, £k,o is the solar luminosity in the K band, 
and the J and K s magnitudes are from the Two Micron 
All Sky Survey (2MASS; mainly the Extended Source 
Catalog, but three galaxies are in the Large Galaxy Atlas 
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of Uarrett et all l2003h and are listed in Table H The 
original relationship uses Johnson K magnitudes, but use 
of K s magnitudes does not introduce significant error to 
the mass. This relation is a linear combination of the T * 
- ( V — J) and - (V — K ) relations presented in Table 1 
of lBell fe de Jond (j2001l) and we subtracted 0.15 dex to 
convert from the scale d Salpeter IMF that the authors 
use to a IKroupal (|200lD IMF. For three galaxies that do 
not have 2MASS data, we used: 

log if = 0.21(5 - K s ) - 1.11 = 0.21(5 - [4.5]) - 1.23, 

(6) 

where we used B-band magnitudes from the Third 
Reference Catalogue of Bright Galaxies (RC3; 
Ide Vaucouleurs et al.l 119911 also list ed in Table Qjj . 
This equation is s imilar to that used inlLee et al.l (2006), 
but we used the iBell fc de Jo"ngl (|2001h T a ble 1 model 
results, we have converted to a IKroupal ()200lD IMF, 
and we assumed (K s — [4.5]) = 0.58. This color is the 
average value from the eleven galaxies in our sample 
where the 2MASS aperture radius ("r_m_ext") and our 
aperture SMA (£> 2 s/2) differ by less than 40% of our 
SMA. Six galaxies have spuriously large/red (K s — [4.5]) 
colors because the 2MASS aperture radius is much 
smaller than our aperture radius (as small as 20% of 
our a perture radius; for a similar effect, see IDale et all 

We then calculated the stellar mass within the 21"- 
diameter circular aperture and the total stellar mass 
with: 

log (M./M ) = logT?+log(L4.B/i4.6,©)-0.4(X.-[4.5]), 

(7) 

as in iLee et all (|2006f ) . To calculate the luminosity of 
the galaxy, L4.5, we used either the 4.5 /zm flux density 
within the 21"-diameter circular aperture or the total 
4.5 /xm flux density from Table [2] the distances presented 
in Tabic [TJ th e absolute 4.5 /mi magnitude of the Sun, 
M 4 . 5 ,0 = 3.3 (jLee et alj|200a lOh et al.ll2008D . and the 
zero-mag nitude flux density of the 4.5 /im band, F° 5 = 
179.7 Jy (jReach et al.ll2005f ). This equation uses the fact 
that (K — [4.5]) ~ for the Sun. We again assumed 
(K s — [4.5]) = 0.58, for all galaxies in this case. For the 
stellar mass surface density (£*), we then divided by the 
deprojected area of the 21"-diameter circular aperture. 

The average offset between masses computed with 
Equations [5] and IH] is 0.2 dex, with masses computed with 
Equation [6] generally smaller. This offset is smaller than 
our estimate of the uncertainty in the mass, described 
further below. The masses computed with Equation [7] 
above are in good agre ement with masses computed us- 
ing the lOh et al.l (|2008l) relationship between K-band and 
4.5 /im band mass-to-light ratios, derived from stellar 
p opulation synthesis mo deling . 

IBell fc de Jonl (|2001f ) cite uncertainties of 0.1 -0.2 dex 
in their mass-to-light ratios, which includes model and 
dust uncertainties and allows for small star formation 
bursts. The uncertainty introduced to the stellar mass by 
uncertainties in L4.5 is of similar order, mainly due to dis- 
tance uncertainties. We assigned a general uncertainty 
of 0.3 dex to the stellar masses, which is the quadrature 
sum of t he above unce r tainti es. This is in agreement 
with the IConrov et al.1 (|2009t ) estimate of typical stel- 
lar mass uncertainties from stellar population synthe- 
sis modeling. The uncertainty in the stellar mass sur- 



face density is dominated by the contribution from the 
mass-to-light ratio. We therefore assigned a typical un- 
certainty of 0.2 dex. 

We estimated a representative oxygen abundance for 
each galaxy to help assess the accuracy and precision of 
our molecular hydrogen and SFR surface densities (Sec- 
tions 13.11 and 13.21) and a lso so we can compare our data 
to the Krumholz et"aT] (|2009l) model (Section O]). We 
estimated the oxygen abundance with the total stellar 
masses calculated from the 4.5 /im flux density a nd the 
mass-metallicity relation of iTremonti et ah| (|2004f) (their 
Equation 3), which assumes a IKroupal (|200lfi " IMF, con- 
sistent with the rest of our analysis. We derive an average 
uncertainty in 12 + \og(0/H) of 0.11 dex. 

3.4. Stability Parameters and Mid-Plane Pressure 

In Section l4~3l we study trends between SFE and sta- 
bilit y, parametrized by a genera lized Toomre Q parame- 
ter iToomrdllQel lRafikovll200lD . We calculated the sta- 
bility parameters within the 21"-diameter circular aper- 
ture. The stability parameter for the gas and stellar com- 
ponents of the disk are 



Qgas — 



ttGS 



, and 



gas 

7rGS* ' 



(8) 
(9) 



respectively, where k is the epicyclic frequency and a gas 
and <7*. r are the gas and radial stellar velocity dispersions. 
iRafikovl ((2001) calculated the stability parameter that 
includes both gas and stars in a thin rotating disk: 



R- 



Qgas+stars Qstars 1 H~ ^ Qgas \ -\- C[ R 



(10) 



where q = ka^^/n, R — cr gas /c*,ri and A; is a free 
parameter that represents the wavenumber of the per- 
turbation. In all cases, the instability condition is 
when Q < 1 (except, str ictly speaking, Q sta rs < 1-07). 
iRomeo fc Wiegertl (|2011l ) included disk thickness in the 
stability estimation by incorporating a factor related to 
the ratio of the ve rtical to r adial velocity dispersion, but 
we settled on the Rafikov (2001) parameter for ease of 
comparison with other works. 

We generally made simila r assu mptions for the com- 
ponents of Q as lLerov et al.l (|2008l ). We assumed cr gas = 
llkms -1 , which is appropriate for the warm neutral 
medium, <7*, r = 1.67er* :Z , and cr*^ = (2ir G U £*/7.3) - 5 , 
where <j*. z is the vertical stellar vel ocity dispersion 
and I* is the stellar scale length (iTamburro et al ' 
2009t iShapiro et al l 120031: Ivan der Kruit fc Searld Il98ll 



van der Kruitl ll988). This cr* equation assumes that the 



disk is isothermal in the z direction, is constant as 
a function of radi us, and h/h* — 7.3 (as observed by 
iKregel et al.l 12002). where ft.* is the stellar scale height. 
Note that the latt er two assumptions are not in conflict 
with a sample like lDalcanton et al.l (|2004[) because those 
authors found only a transition in dust scale height with 
circular velocity, not in stellar scale height. We used the 
scale lengths and stellar mass surface densities from Sec- 
tions [2A4] and [3?3] to derive er* ir values between 9 and 
78kms" 1 . 

For Egas, we multiplied the total- hydrogen surface den- 
sity (£hi+h 2 ) °y a factor of 1.36 to include helium. For 
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Fig. 3. — SFR surface density versus atomic hydrogen surface 
density. Red squares represent galaxies with i> c j rc < 120 kms . 
Blue triangles represent galaxies with v c i IC > 120 km s -1 . The 
circular velocities are from H I rotation curve fits (Paper I). A 
representative error bar is shown in the lower right corner (see 
Sections 13.11 and 13.21 for details). The vertical dashed line is at 
9 Mq pc~ 2 , the typical maximum density for atomic hydrogen. For 
comparison, the s mall black points represent measurements from 
Bigicl ct al. (2008) from seven spiral galaxies sampled at 750 pc 
resolution. 

£», we used the value derived in Section [3.31 For the 
epicyclic frequency, we used the value from Section l2.2.2l 
evaluated at 5.25". For the wavenumber of the perturba- 
tion, k = 2tt/X, where A is the wavelength of the pertur- 
bation, we used the common method of varying A to find 
the minimum (X aa + S t ara (Q ga . s + s t, fl r S .min ) for the region 
(e.g.. lYang etafll2007t lYim et al.ll2011l) . This typically 
results in smaller, less stable <3 ga s+stars compared to Q gas 
and Qstars- We found A = 0.5 - 3.8 kpc at Q gas +stars,min- 
We propagated the uncertainties in Av, t*, Sg as , and 
to derive a typical uncertainty in Q gas of 30%, in Q s tars 
of 60%, and in Q gas+Stars of 40%. 

We also calculated the mid-plane pressure (Ph) with 
the following equation from lElmegreen] (|1989l) : 



Ph ~ ~z GSg as (Sg as H - — S*) 



(11) 



We propagated the uncertainties in l t , S gas , and £* to 
derive a typical uncertainty in P h of 40%. 

4. RESULTS 

4.1. Bulgeless Disk Galaxies on the Kennicutt- Schmidt 

Law 

In this section, we determine whether our galaxy sam- 
ple follows the various versions of the Kennicutt-Schmidt 
law, i.e. the relation between the surface density of gas 
(atomic, molecular, or the sum of both) and SFR (Sec- 
tions 14.1.11 - 14.1.31) . In Section 14.21 we use these results 
to show that there is no transition in SFE at a circular 
velocity of 120 km s" 1 , or at any other circular velocity 
probed by our sample (u c irc = 46 — 190 kms -1 ). The 



Fig. 4. — SFR surface density versus molecular hydrogen surface 
density. Symbols are as in Figure [3] The solid line shows the 
IBieiel et al. (2008) fit to the molecular hydrogen star formation 
law: SgFR °c Sh , . W e discuss the offset of our data relative to 
this fit in Section 14.1.21 



circular velocities were derived from H I rotation curve 
fits in Paper I and we include them in Table [T] for conve- 
nience. 

4.1.1. Atomic Hydrogen Kennicutt-Schmidt Law 

Figure[3]shows the relationship between the star forma- 
tion rate surface density and the atomic hydrogen surface 
density within the 21"-diameter circular aperture, which 
corresponds to physical diameters of 0.7 - 3.2 kpc. Red 
squares denote galaxies with w C j rc < 120 km s" 1 and blue 
triangles denote galaxies with ti c ; rc > 120 kms" 1 . The 
vertical dashed line is at 9 Mq pc -2 and represents the 
typical maximu m H I surface den sity observed in most 
nearby galaxies (|Bigiel et al.ll2008l) . For comparison, we 
have shown as small dots the 750 pc-diam eter regions 
from the s even spiral galaxi es studied in Bigi el et all 
(|2008f) and iLerov et all (120081) Consistent with our as- 
sumptions, iBigiel et al.l (|2008|) used the same Kroupa- 
type IMF and do not include He in their gas surface 
densiti es. An important difference between our datasets 
is that lBigiel et all (|2008|) derived their SFR surface den- 
sities from a combination of FUV (1350 - 1750 A) and 
24 /im data while we use Ha and PAH emission. 

The Spearman rank correlation coefficient for our data 
is 0.5. For comparison, the coefficient for the lBigiel et al.l 
(120081) data i s 0.4. Our data are more correlated than the 
IBigiel et"ail (|2008l ) data, primarily because our sample 
includes several galaxies with large H I surface densities. 

4.1.2. Molecular Hydrogen Kennicutt-Schmidt Law 

Figure 2] shows the relationship between the star for- 
mation rate surface density and the molecular hydrogen 
surface dens ity. The symbols are as in Figure [3] We 
convert the IBigiel et~al1 (|2008l) Sh 2 values to Xco — 
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2.8 x 10 20 cm~ 2 (Kkms " 1 ) 1 to mat c h our assumptions. 
The solid line shows the lBigiel et al.l (|2008|) fit, also con- 
verted to the above Xco- The Spearman rank correla- 
tion coefficient for our data is 0.7, incl uding upper limits . 
For comparison, the coefficient for the lBigiel et Pa l. (2008) 
data is 0.8. 

Our sample of bulgeless disk galaxies, both the low- 
and high-w c i rc objects, appears to lie offset to higher 
star format i on ra te surface densities relative to the 
iBigiel et all ((2008) fit. In what follows, we investigate 
the significance of and possible reasons for this offset. 

Under our assumptions (He not included in the gas 
surface densities, Xqo = 2.8 x 10 20 cm~ 2 [K kms - 1 ] -1 , 
and the Kroupa-type IMF), the lBigiel et al.l (|2008ft fit is 
log(E SF R [Mgyr-ikpc" 2 ]) = a + &log(E H2 [M p C - 2 ]) 
with a = -3.3 ± 0.2 and b = 1.0 ± 0.2. We fit for the 
intercept of our data assuming the Bigiel slope of b = 1.0 
and including our upper limits in the fit and find a — 
—3.0 with an rms of 0.3 dex. 

We determine the significance of the offset between our 
sampl e of bulgeless disk galaxies and the [Bigiel et al.l 
(2008) fit to the molecular hydrogen star formation 
law by rand o mly s electing nineteen measurements from 
IBigiel et~ail (120081). allowing multiple selections of the 
same point ([Press et al.l 11992ft . We assume a slope of 
1.0, calculate the intercept, and repeat this process 10 6 
times. We find that the probability of measuring a inter- 
cept greater than or equal to —3.0 is 4 x 10 ~ 6 . 

We can exclude two possible reasons for the offset. 
First, the offset is not likely due to the measurements 
being centr al values because we confirmed that the cen- 
ters of the IBigiel et al.l (|2008l ) galaxies are not offset 
from the general trend (this can also be seen in Fig- 
ure 10 of Bigiel et al. 2008). Second, our assump- 
tion of a single CO-to-H2 conversion factor probably 
does not lead to underestimated molecular surface den- 
sities because our sample includes only one galaxy for 
which Xqq may be underestimated because of a low 
oxygen abundance, as discussed in Section 13.11 Fur- 
thermore, we do not find that galaxies with lower stel- 
lar mass, and by implicati o n low er- (O /H), are more 
offse t from the IBigiel et all (|2008T ) fit. Note, however, 
that ISchruba et al.1 (|2011f ) did observe that lower-oxygen 
abundance (down to 12 + \og{0/H) = 8.25) galaxies are 
offset to higher molecular SFE in their study of 33 galax- 
ies with directly-measured oxygen abundances. They dis- 
cussed that more observations are needed to determine if 
the offset is due to Xco variation or if it represents true 
SFE variation. 

The quoted IBigiel et al.l (|2008h intercept uncertainty 
(0.2 dex) takes into account variations in star formation 
tracers, uncertainty introduced by estimating the CO(l- 
0) line intensity from CO(2-l) data, and scatter in the 
data. The first of these is particularly relevant in compar- 
ing our datasets be cause we trac e star formation with Ha 
and PAH data and IBigiel et al.l (|2008l ) trace star forma- 
tion with FUV and 2 4 p d ata. Our intercept is larger 
than the Big iel et all (|2008l ) intercept by only 1.2er, if 
we take their error bar as a. In summary, ou r data are 
significantly offset from the IBigiel et al.l ((2008) fit to the 
molecular hydrogen star formation law in terms of statis- 
tical uncertainties, but are nearly consistent if systematic 
e rrors are included . 

iKennicutt et al.1 ([2007) discussed that offsets are ex- 



pected between star formation laws derived from ob- 
servations at different spatial resolution if the power- 
law index of the star formation law is not equal to one 
(Esfr oc ^gas with N ^ 1). In transitioning from high 
to lower resolution, the SFR and gas surface density 
will be decreased by approximately the same factor. If 
N = 1, the lower resolution data will still lie on the 
same relation as the higher-resolution data, but at lower 
gas and SFR surface densities. In contrast, if N > 1 
(N < 1), the lower resolution observations will be pos- 
itively (negatively) offset from the higher resolution ob- 
servations. This effect could contribute to our observed 
offset if TV > 1 because our dat a probe up to 3.2 kpc 
scales and the Big iel et ah! (|2008|) data has 750 pc reso- 
lution. We cannot provide a firm explanat i on for the off- 
set between our data and the lBigiel et al.l (j2008h fit, but 
possible reasons for the offset include the use of different 
star formation tracers and resolutions, Xco variation, 
and true SFE differences. 

4.1.3. Total Hydrogen Kennicutt- Schmidt Law 

Figure[5]shows the relationship between the star forma- 
tion rate surface density and the total hydrogen surface 
density. The vertical dashed line again shows the typical 
maximum H I density that is observed in most nearby 
galaxies. The Spear man rank correlati on coefficient for 
our data and also the lBigiel et al.l (|2008h data is 0.8. This 
correlation is stronger than the correlation between SFR 
surface density and either atomic or molecular hydrogen 
surface density. In their st udy of star formation in the 
atomic-dominated regime, ISchruba et all (|20T1 found 
that the rank correlation coefficient between SFR sur- 
face density and total hydrogen surface density is similar 
to or somewhat larger than the value for the correlation 
with the molecular hydrogen surface density. However, 
they noted that data can be correlated with a strong 
rank statistic even if the parameters are related by dif- 
ferent functions in the atomic- and molecular-dominated 
regimes. As in the Schruba et al. study, our stronger 
total hydrogen correlation is not likely related to funda- 
mental physics. 

The solid black line shows the iKrumholz et al.l (|2009h 
model with 0.4 Zq metallicity. The model assumes a 
clumping factor, c, that is the inverse of the filling fac- 
tor of ~ 100 pc-sized atomic-molecular complexes in the 
beam. This value is not constrained well by dat a, so we 
assume the same value as IKrumholz et aT . (2009): c = 5 . 
The magenta dot-dashed line shows the lKennicuttl ()1998[ ) 
fit to a sample of normal and starburst galaxies, where 
the measurements are averages over t he entire optical 
disk. The blue dotted line shows the IKennicutt et al.1 
(2007) fit to regions in M51 t hat w ere studied at 520 pc 
resolution. IKennicutt et al.l (|2007l ) attribute the offset 
between the blue and magenta lines to dilution effects 
when the power-law index of the star formation law is 
not 1.0 (as discussed in Section H.1.2p . 

4.2. No Transition in SFE with Circular Velocity 

In Section [TJ we discussed that there could be a SFE 
transition in bulgeless disk galaxies at v c j rc = 120 km s -1 , 
depending on the star formation model assumed. In this 
section we show that there is no offset between the low- 
and high-u c ; rc galaxies on the star formation law. In Sec- 
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F ig. 5. — SFR sur face density versus total hydrogen surface density. Symbols are as in Figure [3] The magenta dot-dashed line shows 
the Kcnnicutt (1998) fit to a sa mple of 61 normal spir als and 36 star burst galaxies, where the measurements are averages over the optical 
disk. The dotted line shows the [Kcnnicutt ct al. (2007) fit to M51 data at 520 pc resolution. The solid black line shows the Krumholz ct al. 
(2009) model at 0.4 Z Q metallicity. 



tionO we interpret this result to discuss the relationship 
between SFE and the scale height of the cold ISM. 

In Figures [3l |U and [5l there is a slight trend for low- 
■fcirc galaxies to lie at lower surface densities compared 
to high-u c ; rc galaxies, with average low- versus high-z; c i rc 
surface densities as follows: log(SHi) = 0.9 versus 1.0, 
log(EH 2 ) < 0.6 versus = 0.9, log(£Hi+H 2 ) = 1-1 versus 
1.2, and log(SsFR) = —2.4 versus -2.2. However, the 
differences are not significant; the Kolmogorov-Smirnov 
(K-S) test probabilities that the low- versus high-w c i rc 
galaxies are drawn from the same population in £r- 2 , 
£hi+h 2 , an d Ssfr are 0.04, 0.2, and 0.2, respectively. 

Figure [5] shows the distribution of SFE in the low-w circ 
(red solid) and high-i> c i rc (blue dashed) galaxies. If the 
slope of the star formation law is not 1.0, a sample that 
follows the relation would have SFE that varies as a func- 
tion of Shi-(-h 2 , which will tend to broaden the distri- 
butions. However, our low- and high-w c i rc objects have 
similar £hi+h 2 distributions, so this should not lead to 
spurious offsets between the samples. Furthermore, the 
£hi+h 2 range is not so great that it would hide any dif- 
ferences between the samples. We find no significant dif- 
ference between the SFE distributions, with a K-S test 
probability of 0.8 that the low- and high-v c i rc galaxies 
are drawn from the same population. 

We conclude that our sample of bulgeless disk galax- 
ies does not show a strong transition in SFE at w C i rc = 




-10 -9.8 -9.6 -9.4 -9.2 

log(SFE) = log(2 srB /S 

HI + H„ [yr-D 

Fig. 6. — SFE histograms of the low-n c ; rc (red solid) and high- 
v c i IC (blue dashed) samples. There is no significant difference be- 
tween the distributions. 
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Fig. 7. — SFE versus circular velocity. There is no circular 
velocity at which a transition in SFE occurs. 

120 km s -1 , where iDalcanton et al.l (|2004f ) concluded 
that there is a strong transition in dust scale height. In 
Figure [7] we show that there is no circular velocity at 
which there is a transition in SFE. Furthermore, there 
is no strong trend of SFE with circular velocity within 
the range we probe. Finally, we find no difference in the 
ratio of the molecular to atomic surface density between 
the low- and high-w circ galaxies (the K-S test probability 
that the samples are drawn from the same i? mo i popula- 
tion is 0.4), nor do we find a transition in i? mo i at any 
circular velocity (Figure |8]). 

4.3. SFE Trends with Stability and Mid-plane Pressure 

In this section, we address whether our sample shows 
a transition or trends in stability or mid-plane pressure, 
using the Q gas , <3 sta rs, Qgas+ S tara,min, and P h values cal- 
culated in Section 13.41 Figure [9| shows the SFE versus 
<3 g as (left), Q stars (middle), and <3 gas +star S: min (right). 
We find no unstable regions. The SFE generally de- 
creases with larger, more stable Q values, although the 
correlation is not strong. For SFE versus <5 ga s, Qstars, 
and Q gas +stars,min, we find Spearman rank correlation 
coefficients of -0.2, -0.3, and -0.2, respectively. The sign 
of this correlation is as expected, but one might instead 
have expected a sharp decrease in SFE when Q rises 
ab ove 1. Thes e resul ts are qualitatively similar to those 
in iLerov et al.l (120081) . 

IDalcanton et al.1 |2004l ) concluded that galaxies with 
Vdrc < 120kms _1 are generally stable while galaxies 
with v c ; rc > 120 km s -1 are generally unstable, especially 
in the central (r < Z„) regions. Contrary to these re- 
sults, we see no evidence for a transition in stability at 
Vcirc = 120kms _1 ; the K-S test probability that the low- 
and high-u c j rc galaxies are drawn from the same popu- 
lation of <5 gas , Q s t a rs, and Qgas+stars.min is 0.5, 0.2, and 
0.8, respectively. 

There are two principal differences between our 



Fig. 8. — Ratio of the molecular to atomic gas surface density 
(iJ mo l) versus circular velocity. There is no circular velocity at 
which a transition in i? mo i occurs. 



assumptions for the stability inputs and those in 
IDalcanton et al.l (|2004D . First, we use a constant gas ve- 
locity dispersion and they used the quadrature sum of the 
velocity dispersion of atomic and molecular gas (10 and 
5kms _1 , respectively), weighted by the relative mass 
surface densities of the components. We use these as- 
sumptions and still find no difference in stability between 
the low - and high-i; c j rc galaxies. Second. IDalcanton et all 
(2004) estimated the molecular hydrogen surface densi- 
ties for their sample from a scaling with circular veloc- 
ity (Eh 2 = (wcirc/47.1 kms -1 ) 2 49 ). This scaling relation 
was derived from a similar sample of galaxies with Eh 2 
values from lRownd &: Youngl (| 19991 ) and w c irc estimated 
from single-dish H I data. Molecular hydrogen surface 
densities calculated using the scaling with w C j rc are too 
large by a factor of six for our sample of high-u c j rc galax- 
ies and too large by a factor of two for our sample of 
low-u c i rc galaxies. This discrepancy leads to smaller, less 
stable <3 g as and Q gas +stars in the high-w c ; rc galaxies rela- 
tive to the low-u c i rc galaxies. However, even under these 
assumptions there is no significant difference between the 
low- and high-u c i rc galaxies in their Q values. Note that 
the difference in our molecular hydrogen surface densi- 
ties compared to the values derived from the scaling with 
u c i rc may also be related to resolution differences. Our 
measuremen t s are within the central 21" while the 
iRownd &: Youn g (1999) data have 45" resolution and a 
similar distance distribution. 

We test whether a stability transition occurs at any cir- 
cular velocity in Figure [101 where we plot <3 gas +stars,mm 
versus t> C i rc . There is no circular velocity below which 
the regions are stable and above which the regions are 
unstable. If any trend is present, it is that higher circular 
velocity objects are more stable. 

We plot the SF E versus the mid-p lane pressure in Fig- 
ure [TlJ As seen in lLerov et al.l (|2008l) . we find an increase 
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Fig. 9. — SFE (Esfr/£hi+h 2 ) versus stability parameters Q gas (left), Qstars (middle), and Q gas +stars, min (right)- The symbols are as 
in Figure [3] Note that the condition for instability is Q < 1, so none of our galaxies are strictly unstable. We describe our derivation of 
the errorbar in Sections 13.21 and 13.41 
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Fig. 10. — Stability parameter, including the gas and stellar 
contribution, versus circular velocity. There is no circular velocity 
at which a transition in stability o ccurs . We describe our derivation 
of the vertical errorbar in Section 13. II 

in SFE with increasing mid-plane pressure, but we do not 
probe high enough P^ values or have enough data points 
to sample the constan t -SFE region of the diagram that 
is clear in iLerov et al.l (|2008| ). We find no difference be- 
tween the low- and high-i> c i rc galaxies in their Ph distri- 
butions, with a K-S test probability of 0.7 that they are 
drawn from the same population. Even focusing only on 
galaxies with the same Shi+h 2 (within the uncertainty 
in the parameter), the median pressure is the same in 
low- and high-u c ; rc galaxies. 



Fig. 11. — SFE versus mid-plane pressure. Symbols are as in 
Figure [3] 

4.4. Dependence of Star Formation on Metallicity 

In this section, we use the oxygen abundance esti- 
mated from the stellar mass to study how star forma- 
tion depends on metallicity and compare the s e resu lts to 
recent theoretical work by Krumho lz et al.l (12009T). At 
a give n total gas surface density, the iKrumholz et al.l 
( 2009) model predicts lower SFE at lower metallicity be- 
cause H2 survival requires a higher column density as 
H2 self-shielding becomes more important than shield- 
ing by dust. This metallicity dependence distinguishes 
the model from other leading models, such as the model 
where mid-plane pressure determines where the ISM is 
molecular. 

The left panel of Figure [12] shows the total hydrogen 
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Fig. 12. — Left: SFR surface density versus total hydrogen surface density, as in Figure[3] but showing a larger range to demonstrate 
the model trends. The data are divided by oxygen abundance, where orange pentagons represent galaxies with 12 + log (O/H) < 8.7 and 
purple circles represent galaxies with 12 + log (O/H) > 8.7. The lines show the Krumholz ct al. (2009) model with a metallicity of 0.1 Zq 
(orange solid) and 1.4 Zq (pur ple dashed), which corr espond to the lowest and highest metallicities in our sample. Right: Histograms of 
the offset of the data from the Krumholz et al. (2009) model with a metallicity of 0.4 Zq. The orange line represents the distribution of 
galaxies with 12 + log (O / H) < 8.7 and the purple line represents the distribution of galaxies with 12 + log (O/H) > 8.7. There is no 
significant difference between the distributions. 



star formation law, as plotted in Figure [5l except the 
galaxies are divided into low-(0/H) and high-(0/ H), 
with the division at 12 + \og(0/H) = 8.7. The or- 
ange solid and p urple dashed lines show the model of 
IKrumholz et~all (|2009D at 0.1 Zq and 1.4 Zq, respec- 
tively. These metallicities correspond to the extrema 
of our dataset, assumi ng Z/Zq = (Q/H)/(Q/H)q an d 
12+log(O/fl") = 8.86 (jDelahave fc Pmsonneaul11l2006l ). 
We see no evidence that the low- and high-(0/if) galax- 
ies cluster towards the low- and high-Z models, respec- 
tively. 

The right panel of Fi gure IT2l shows t h e offse t of our 
data from the 0.4 Z Q IKrumholz etail (|2009D model, 
which provides the best fit to our dataset as a whole 
(although note that we would expect the best fit to be 
the 0.7 Zq model, as that corresponds to the average 
(O/H) of our full sample). The orange solid and purple 
dashed lines show the distributions of our low- (O/H) and 
high-(0/iJ) galaxies, as divided in the lef t pane l. If our 
data clearly followed the IKrumholz et al.1 ()2009f ) model, 
high- (O/H) galaxies would be positively offset from the 
0.4 Zq model and \aw-(0/H) galaxies would be nega- 
tively offset. We find no significant difference between 
the offset of the low- versus high- (O/H) galaxies (the K- 
S test probability that the samples are drawn from the 
same offset population is 0.14). A second order effect 
is that the offset from the 0.4 Zq model is expected to 
decrease with increasing £hi+h 2 • This may be related 
to the narrower offset distribution for the high-(0/ H) 
galaxies, which tend to have higher Ehi+h 2 • 

Figure [T3] shows that there is a correlation between 
SFE and oxygen abunda nce in our data. For com parison, 
the line represents the IKrumholz et al.1 ()2009l ) model, 
where the SFE depends on metal licity, S gas , and the fill- 
ing factor. We have plotted the IKrumholz et al.l (|2009l ) 
model with log(£Hi+H 2 ) = 0.8 (no He). This is a low 
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Fig. 13. — SFE versus oxygen abundance, scaled to the solar 
value. The solid line demonstr ates the SFE trend with metallicity 
for the IKrumholz et al.l d2009l) model at log(£ HI+H2 ) = 0.8. 



Effl+H 2 va hie compared to our data. Higher Ehi+h 2 
models reach constant SFE at lower metallicity. The 
mismatch between the average Ehi+h 2 in the data and 
the best-fit Shi+h 2 used in the model is likely due to our 
assumptions for Xqo and the IMF in the data, our fill- 
ing factor assumption in the model, and to uncertainties 
in the metallicity scale. For the latter, various strong- 
line metallicity methods return 12 + log(0/_ff) values 
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as disparate as 0.5 dex (|Kewlev fc Ellisonl 120081) . and it 
is unknown which calibration best aligns with the solar 
value. Therefore, the model line in Figure [T3] only rep- 
resents the expected trend of lower SFE at lower metal- 
licity, at a given S gas . At this level, the data do show 
the expected trend, with a Spearman correlation coeffi- 
cient of 0.6. However, upon closer inspection, this agree- 
ment is mainly due to the fact that lower- (O/H) galaxies 
tend to have lower Ehi+h 2 (this is most obvious in Fig- 
ure [T^l although also note that the K-S test probability 
is 0.05 that the Ehi+h 2 distributions of the low- and 
high- (O/H) samples are drawn from the same popula- 
tion, which is not strong evidence for a difference). Be- 
cause of their lower Ehi+h 2 values, the \ow-(0/H) galax- 
ies are closer to the atomic-dominated regime of the star 
formation law wher e SFEs are lower . This is of course 
consistent with the iKrumholz et al.l (|2009() model, but 
is not a discriminating test of the model because many 
other propert ies correlate with gas surface density (for 
examples, see lLerov et a l. 2008). A sample with a range 
of directly-measured metallicities within a small range of 
gas surface density would provide a more discriminating 
test . In summary, we see no clear evidence to support 
the IKrumholz et al.l (|2009t ) model, but our data are also 
not inconsistent with it. 

5. DISCUSSION 

In Section 14.21 we found that there is no transition in 
molecular fraction or SFE at any circular velocity probed 
by our sample. In Section 14. 3[ we found that all our 
galaxies are formally stable. While we did find a general 
trend of decreasing SFE with larger, more stable Q val- 
ues, we found no sharp transition in stability at any cir- 
cular velocity. Finally in Section POI we found that SFE 
decreases at lower oxygen abundance, bu t the trend is not 
a part icularly constraining test of the IKrumholz et al.1 
(2009) model. In this section, we first address our as- 
sumption that there is a transition in dust scale height 
at i> c i rc = 120 km s -1 , rather than a transition in dust 
content. We then interpret our results to discuss the re- 
lationship between SFE and the scale height of the cold 
ISM. Finally, we comment on the scale of physical pro- 
cesses that affect star formation and discuss our results 
in the context of leading star formation models. 

5.1. A Transition in Dust Scale Height versus Dust 
Content 

In this section, we investigate the ar gument that the 
dust s tructure transition observed by iDalcanton et al.l 
(|2004l) is due to a transition in dust scale height rather 
than due to a transitio n in the amount of dust present. 
IDalcanton et al.l (|2004|) came to this conclusion because 
the dust structure transition occurs over a narrow range 
in circular velocity, where a large change in the DGR, 
and therefore dust content, is unexpected. All but one 
galaxy in our sample are detected at Infrared Astronomi- 
cal Satellite 60 and 100 /im, which indicates that there is 
at least some dust in the low-w c i rc galaxies. We estimated 
a DGR proxy as the ratio of the total infrared lumi- 
nosity ^Lrii^^cjdj^lated with IR AS 25, 60, and 1 00 /xm 
data (jMoshir et al.lll990Q and the iDale fc Heloul (|2002l ) 
relation, to the combined atomic and molecular hydro- 
gen mass. We find no transition in this DGR proxy at 
^circ = 120kms _1 , nor do we find any correlation be- 



tween the DGR proxy and circular velocity. Ltir is not 
necessarily proportional to the dust mass because the 
dust temperature may not be the same for all the galax- 
ies; nevertheless, there is clearly a significant amount of 
dust in the low-w c i rc galaxies. 

Two recent studies have found convincing evidence 
that low-y c irc sp irals have large dust scale heights. 
ISeth et all ([2005D carried out a resolved stellar popula- 
tion study of six edge-on, late-type spirals with v c i IC = 
67 — 131kms _1 and found that the scale height of the 
young (< 10 8 yr) stellar population, which presumably 
form ed from the cold ISM, is larger than in the Milky 
Way. iMacLachlan et all (|2011l ) modeled the spectral en- 
ergy distribution (SED) of three edge-on, low surface 
brightness galaxies with v c j rc = 88 — 105 km s^ 1 and 
found that a significant amount of dust must be present 
to account for the FIR (70 and 160 /im) emission, but 
the dust must have a large scale height such that it does 
not significantly obscure the optical emission. The au- 
thors concluded that the galaxies have dust scale heights 
greater than or equal to the stellar scale heights. This 
is in contrast to modeling studies of high surface bright- 
ness galaxies, which concluded that the du st scale height 
is ab out half the stellar scale height (e.g.. Kilouris et ah! 
fl999h. 

iHunter fe Elmegreenl (|2006l ) provided an alternative 
ex planation for the dust structure transition observed 
by IDalcanton et all (|2004[) . The authors studied a large 
sample of irregular galaxies and found that the average 
B-band surface brightness is smaller than in higher-mass 
spiral galaxies. Based on this result, they suggested that 
lower stellar surface density is the cause of the larger 
scale heights in low-w circ galaxies. They discussed that 
disk stability could correlate with the dust structure but 
not be the cause of the transition because the gas scale 

height (/ig as oc -g g-, where p is the m ass volume density 
of gas and stars; iBlitz fc Rosolowskvll2004l ) shares many 
of the same parameters with the stability parameter. 

The authors discussed that dust opacity may also con- 
tribute to the observed difference in dust structure be- 
cause of two effects. First, lower-w C i rc objects should have 
lower metallicity, and therefore lower dust content. Sec- 
ond, the scale length of a lower-mass galaxy is smaller 
and therefore a low-v c i rc , edge-on galaxy will have less 
depth from the edge to the center over which to accumu- 
late dust column density compared to a high-i> C i rc galaxy. 
Note that no scale height transition is needed in this in- 
terpretation. However, dust opacity is not likely the sole 
cause of the dust structure transition because of the rea- 
sons listed above. IHunter fc Elmegreenl ()2006l ) suggested 
that the observed dust structure transition is likely due to 
a combination of stellar surface density and dust opac- 
ity. The authors agreed that a scale height transition 
does occur, so if this interpretation is correct, we can still 
constrain the effect of scale height differences on SFE. 

In Sections [5.21 and |5"31 we assume that the dust scale 
heights in low-z) C j rc galaxies are a factor of t wo la rger 
than in high-w c i rc galaxies. IDalcanton et al.1 (|2004l ) es- 
timated this factor by examining the dust morphology 
in Hubble Space Telescope images of a couple edge-on, 
late- type galaxies. This factor is consistent with the 
SEP modeling resul ts of IMacLachlan et al.1 fj201 lh and 
IXilouris et al.l (|1999f ) if the stellar scale height distribu- 
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tions of lo w- and high-t; c i rc galaxie s have significant over- 
lap, which IDalcanton et al.l ()2004f ) found to be true (al- 
though there are low-i; c i rc galaxies where the stellar scale 
height is about half the value in high-v c j rc galaxies, in 
which case the dust scale heights may be comparable). 

5.2. The Relationship between SFE and Scale Height 

In Section l4.21 we found no transition in SFE at u c ; rc = 
120 km s _1 or at any circular velocity probed by our sam- 
ple. In this section, we interpret this result to discuss the 
relationship between SFE and the scale height of the cold 
ISM. For this discussion, we make a number of assump- 
tions. First, we assume that our sample of bulgeless d isk 
galaxies is similar to that of IDalcanton et al.l (|2004f) in 
that galaxies with u circ > 120 km s^ 1 have narrow dust 
lanes while galaxies with v c j rc < 120 km s~ x have no obvi- 
ous dust lanes. This is a reasonable assumption because 
it was our main consideration in choosing the sample, but 
because our galaxies are moderately inclined rather than 
edge-on, we are unable to d irectly measure the ver tical 
dust structure. Second, as in IDalcanton et all (|2004l) . we 
assume that the dust scale heights in the low-w C i rc galax- 
ies are larger than in the high-u c i rc galaxies (Section 15.11 
addresses this assumption). Finally, we assume that the 
molecular gas and dust s cale heights are comparable. 

IDalcanton et al.l (|2004l ) suggested that there may be 
a SFE transition associated with the dust scale height 
transition at w c ; rc = 120 km s -1 . The authors gave an 
example that assumes the true star formation law is a 
correlation between the volume density of gas (p g as) and 
the SFR volume density (psfr)' Psfr oc p^ s . The larger 
scale heights of low-w circ galaxies lead to lower gas vol- 
ume densities. Depending on the index, TV, a low-w c i rc 
galaxy with the same gas surface density as a high-u c i rc 
galaxy can have a lower SFR surface density and there- 
fore a lower SFE. To consider this point, we first assume 
the same E gas for a low- and high-w circ galaxy and the fol- 
lowing relationships between the surface and volume den- 
sity of gas and SFR: E gas oc p gas h, and E S fr oc p S FR h, 
where h is the scale height of the star forming gas and 
newly formed stars and the proportionality constant is 
the same for both relationships. We set f3 equal to the ra- 
tio of dust scale heights in l ow-t> c i rc (lv) versus high -t> c j rc 
(hv) galaxies: P = hi v /hh v - IDalcanton et al.1 (j2004f ) very 
approximately estimated this ratio to be about two. Un- 
der these assumptions, Esfr.Iw = P N ~ X Esfr,1v, where 
EsFR,hv and Esfr,1v are the star formation rate surface 
densit ies of the low- and high -ti r i rr , galaxy with the same 
E gas . IDalcanton et al.l (|2004[) discussed the case where 
P = 2 and N = 1.5, where we expect the high-u c i rc galaxy 
to have a Esfr that is a factor of 1.4 larger than the 
low-v c i rc galaxy. Note that if N = 1 there is no expected 
difference in star formation rate surface density. 

In the star formation law plots of Figures [3l 21 and O 
we are sensitive to offsets in intercept between the low- 
and high-w c j rc samples that are greater than the uncer- 
tainty in the intercept, which is ~ 0.3 — 0.4 dex, assum- 
ing an uncertainty in Esfr of 0.2 dex and the number 
of galaxies in ou r low- and high-^ r i rr , s amples. In the 
case discussed by IDalcanton et ail ()2004f ). we expect the 
low-w c i rc galaxies to be offset to lower Esfr by 0.15 dex. 
T herefore, we cannot exclu de offsets at the level expected 
bv IDalcanton et al.l (|2004l) . 

The star formation law assumed above is not likely cor- 



rect given that recent studies have found a strong molec- 
ular star formation law and no atomic gas star forma- 
tion law. However, we can also determine the expected 
Esfr offset if the molecular fraction is set by the mid- 
plane pressure and the molecular SFE is constant (see 
also Section [lj. We assume the same Ehi+h 2 f° r a l° w - 
and high-u circ galaxy, i? mo i oc P h , and P h cx ,ohi+h 2 cr gas , 
where phi+h 2 is the total hydrogen volume density and 
the velocity dispersion of the gas is constant. In this 
scenario, /9hi+h 2 , Ph, and i? mo i are lower by a factor 
of P in the low-w c i rc galaxy compared to the high-ii c i rc 
galaxy. The expected offset in Esfr depends on i? mo i, 
which varies from 0.1 to 2.7 in our sample. If ft = 2 
and i?moi of the high-w c i rc galaxy is 0.1 (2.7), we expect 
Esfr to be lower by 0.3 dex (0.1 dex) in the low-u c i rc 
galaxy compared to the high-w C i rc galaxy. Given our un- 
certainties, this level of offset would also be difficult to 
detect. However, we can reject our assumption that Ph 
is lower by a factor of two in low-i> c i rc galaxies. This is 
evident from Section 14. 3[ where we found no significant 
difference in the mid-plane pressure distributions of the 
low- and high-w C i rc galaxies, even when considering only 
objects with the same Ehi+h 2 - Furthermore, we found 
no difference in the i? mo i distributions of the low- versus 
high-w circ galaxies (Section r4~2")l . 

The two scenarios explored above predict offsets in in- 
tercept between the low- and high-w C i rc samples that are 
less than the uncertainty. Therefore, we cannot clearly 
exclude these options. Nonetheless, our data show no 
evidence for a strong transition in SFE at any circular 
velocity. A simple interpretation that is consistent with 
our data is that low-i> c i rc galaxies have a lower num- 
ber of molecular clouds per unit volume compared to 
high-?; c i rc galaxies at the same Ehi+h 2 , but lower only 
by the ratio of the cold ISM scale heights in low- ver- 
sus high-u C j rc galaxies. This results in the same total 
number of molecular clouds within the beam for a low- 
and high-u C j rc galaxy at the same Ehi+h 2 and thus gives 
the same Eh 2 and Esfr (assuming the molecular clouds 
have the same density and that we average over evo- 
lutionary effects). Note that the above applies in the 
molecular-dominated regime. We have few data points 
below Ehi+h 2 ~ 9Af Q pc~ 2 , but expect that Esfr can 
vary substantially for galaxies with the same Ehi+h 2 j de- 
pending on the physical processes that affect the molecu- 
lar fraction (e.g., those processes discussed in Krumholz 
et al. 2009 and Ostriker et al. 2010). 

In conclusion, we interpret our result that there is no 
transition in SFE at any circular velocity as evidence that 
scale height differences at the level of about a factor of 
two do not significantly affect the molecular fraction or 
SFE in bulgeless disk galaxies. However, offsets in SFE 
below our uncertainty level are still possible. 

5.3. Comparison to Star Formation Models 

IDalcanton et al.l (|2004f) very approximately estimated 
that the dust scale heights of low-w C j rc galaxies are about 
a factor of two larger than high-u c ; rc galaxies. Assuming 
our sample has a similar range in scale height, our results 
indicate that these scale height differences, which lead to 
gas volume density differences also at the level of a fac- 
tor of about two, do not lead to obvious differences in 
the SFE. Our results favor star formation models where 
small-scale physical processes are more important than 
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processes that act on larger scales, of order the dust and 
cold gas scale height (10s to 100 pc). Based on their com- 
parison to ma ny star formation models without an ob- 
vious favorite . iLerov et al.l (2008) discussed that physics 
below their resolution of 750 pc is likely most important 
for determining the SFE. We contribute with a further 
constraint that the SFE is likely affected primarily by 
processes that act on scales smaller than the cold gas 
and dust scale height. 

We cannot exclude all star formation models that in- 
clude large-scale physics because there are processes that 
may affect star formation but are neither affected by 
nor affect the scale height. For example, our sample 
has no power to constrain the effe cts of large-scale ra - 
dial processes, like shear (see, e.g.. iHunter et alj fl998). 
on star formation. In addition, star formation may be 
affected by environmental properties that depend on the 
gas volume density but also depend on other variables 
that counteract a variable volume density. For exam- 
ple, the pressure in the ISM is related to the gas volume 
density and velocity dispersion: P cx PgasO"g as . While 
we expect the volume density to be lower in low-u c i rc 
galaxies, the velocity dispersion may be larger. With the 
right combination of p gas and a gas , there could be no dif- 
ference in pressure between low- and high-u C j rc galaxies. 
While the latter argument illustrates our limitations, we 
note that there is currently no strong observational rea- 
son to assume different gas veloc ity dispersions b e tween 
the low- and high-« c ; rc galaxies. iTamburro et al.l (2009) 
found some variation in the central H I velocity disper- 
sion in eleven galaxies ranging from early-type spirals to 
irregulars, but more observations are needed. 

Would leading star formation models have predicted a 
difference in SFE in galaxies with scale heights tha t differ 
by about a factor of two? In the lKrumholz et"al"1 (2009) 
model, EgFR is a function of metallicity, S ga s, and the 
beam filling factor of ~ 100 pc-sized atomic-molecular 
complexes. There is no direct dependence on the scale 
height, so we would not expect a transition in SFE unless 
there is a transition in the metallicity or filling factor with 
scale height. We find no transition in oxygen abundance 
at any circular velocity in our data, but there should be a 
correlation between thes e two properties given the mass- 
metallicity relation (e.g.. lTremonti et al.l l2004). which we 
also do not see. The filling factor of star-forming com- 
plexes is not well constrained, although there is no a pri- 
ori reason to suppose th at it would be different in low- 
versus high-w c i rc objects. iKrumholz et al.l (|2009l ) do not 
predict a difference in SFE in galaxies with different scale 
heights and the fact that we did not find a transition in 
SFE at u c hc = 120 km s" 1 is not a strong constraint on 
this model. 

In the model where the mid-plane pressure sets the 
molecular fraction and the molecular SFE is constant, 
the pressure is proportional to the gas volume density, 
which we expect to vary between the low- and high- 
w c irc galaxies. If the gas velocity dispersion is fixed, we 
would expect the low-u c i rc galaxies to have lower molec- 
ular to atomic surface density ratios and lower SFEs rel- 
ative to the high-v c j rc galaxies . However, we found no 
difference between the mid-plane pressure, molecular to 
atomic surface density ratio, or SFE distributions of the 
low- and high-u c ; rc galaxies. If our assumptions are cor- 



rect, our result is inconsistent with this model. How- 
ever, the offset in SFE expected for galaxies with cold 
gas scale heights that differ by about a factor of two may 
be less than our SF E uncertainties. Furthermore, the 
lOstriker et al.l ()2010l ) model relates the molecular frac- 
tion (or in their terms, the fraction of gas in gravita- 
tionally bound complexes) to the pressure of the diffuse 
component of the ISM. We have a constraint only on the 
scale height and volume d ensity of the co l d com ponent 
of the ISM; therefore the lOstriker etail (|2010D model 
may not predict molecular fraction and SFE differences 
in our sample. In general, our results are somewhat more 
consistent with l ocal m odels of star formation, like the 
IKrumholz etHI (|2009T) model, but we do not find con- 
clusive ev idence for or against e ither the IKrumholz et al.l 
(|2009T ) or lOstriker et alJ (|2010t l model. 

One final matter to address is whether central measure- 
ments are sufficient to determine if there is a transition 
in molecular fraction, SFE, and/or stability at the dust 
structure transition of w c ; rc = 120 km s" 1 . One might 
question the use of central measurements because dust 
structure, SFE, and stability may be affected by the 
higher gas and stellar densities and shorter dynamical 
times characteristic of these regions. The only approach 
that will fully address this concern is to obtain off-center 
measurements of the above properties. Our single-beam 
CO data currently limit us from carrying out this anal- 
ysis. Meanwhile, there is some evidence that our central 
pointings are sufficient to address these questions. First, 
our measurements trace a significant fraction of the disk: 
the 21" aperture probes physical scales of 0.7 — 3.2 kpc, 
which is similar to the disk scale lengths of our sample 
(0.7 — 3.4 kpc). Second, the central regions of bulgeless 
galaxies are more morphologically and kinematically sim- 
ilar to the outskirts than in a galaxy with a bulge. Fi- 
nally, we expect a galaxy to be less stable against grav- 
itational collapse in the center compared to the outer 
disk because the gas and stellar surface densities are 
larger. However, we find that both the low- and high- 
fcirc galaxies are stable. This suggests that there would 
also not be a stability transition at u c i rc = 120kms _1 
in off-center measurements because both the low- and 
high-v c h- c galaxies would be more stable. 

6. SUMMARY 

We have presented a study of star formation in twenty 
moderately-inclined, bulgeless disk galaxies. We found 
no transition in star formation efficiency (Esfr/£hi+h 2 ) 
or disk stability at w C i rc = 120 km s -1 . This circular ve- 
locity was previously found to be associated with a tran- 
sition in the vertical dust structure of edge-on, bulge- 
less disk galaxies that is most likely due to a transition 
in the scale height of the cold ISM. We also found no 
transition in star formation efficiency or disk stability at 
any circular velocity probed by our sample. Our results 
demonstrate that the scale height of the cold ISM does 
not play a major role in setting the molecular fraction or 
the star formation efficiency. We also found decreasing 
star formation efficiency with lower oxygen abundance, 
which we estimated from the stellar mass. This result is 
consistent with the recent IKrumholz et al.1 (|2009t ) model, 
but a sample with a large range of metallicities within a 
small range of gas surface density would provide a more 
constraining test of the model. In general, our results 
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are most consistent with local models of star formation 
that include physical processes that act on smaller scales 
than the dust and cold gas scale height (10s to 100 pc). 
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(4) 


(arcscc) 
(5) 


B 
(mag) 
(6) 


Mb 
(mag) 
(7) 


Type 
(8) 


1/ 

(kms -1 ) 
(9) 


(kms -1 ) 
(10) 


PA 

(°) 
(11) 


(°) 
(12) 


Won 

vv 20 
(kms -1 ) 
(13) 


NGC 0337 


00 


59:50.0 


-07:34:41 


20.7 [T88] 


173 


11 


44 


-20.14 


7 





1646 ± 2 


145+4 


118 ± 5 


44 + 2 


261 


PGC 3853 


01 


05:04.8 


-06:12:46 


11.4 [T08] 


250 


11 


98 


-18.30 


7 





1094.7 ± 0.4 


128. ltj 8 


105.3 ± 0.2 


41.4 ± 1.1 


192 


PGC 6667 


01 


49:10.3 


-10:03:45 


24.6 [T88] 


173 


12 


92 


-19.03 


6 


7 


1989.2 ± 0.6 


155±* 


122.9 ± 1.8 


34.0 ± 1.1 


198 


ESO 544-G030 


02 


14:57.2 


-20:12:40 


13.9 [T08] 


123 


13 


25 


-17.47 


7 


7 


1608.4 ± 1.0 


100. 9+. 1 ' 4 


107.6 ± 1.1 


48.5 ± 1.2 


146 


UGC 1862 


02 


24:24.8 


-02:09:41 


22.3 [T08] 


99.6 


13 


47 


-18.27 


7 





1382.9 ± 0.4 




21.7 ± 1.7 


43 + 4 


125 


ESO 418-G008 


03 


31:30.8 


-30:12:46 


23.6 [T08] 


70.5 


13 


65 


-18.21 


8 





1195.4 ± 0.3 


74.i±g;| 


317.9 ± 1.1 


55.6 ± 1.4 


140 


ESO 555-G027 


06 


03:36.6 


-20:39:17 


24.3 [T88] 


138 


13 


18 


-18.75 


7 





1978.7 ± 0.4 


190± 4 o 


221.5 ± 0.3 


21 + 4 


162 


NGC 2805 


09 


20:20.4 


+64:06:12 


28.0 [T88] 


379 


11 


17 


-21.07 


7 





1732.6 ± 0.6 


81±I 


300 ± 3 


38 + 4 


120 


ESO 501-G023 


10 


35:23.6 


-24:45:21 


7.01 [T08] 


208 


12 


86 


-16.37 


8 





1046.8 ± 0.7 




224 ± 2 


37 ± 12 


83 


UGC 6446 


11 


26:40.6 


+53:44:58 


18.0 [T08] 


213 


13 


30 


-17.98 


7 





645.5 ± 0.6 


' y - ' -0.8 


189.4 ± 0.5 


52.5 ± 1.9 


150 


NGC 3794 


11 


40:54.8 


+56:12:10 


19.2 [T08] 


134 


13 


23 


-18.19 


6 


5 


1384.9 ± 0.7 


103.3+J-J 


123.1 ± 1.1 


54.8 ± 1.3 


182 


NGC 3906 


11 


49:40.2 


+48:25:30 


18.3 [...] 


112 


13 


50 


-17.81 


7 





959.44 ± 0.7 


65+ 30 


180 ± 20 


16 + 5 


49 


UGC 6930 


11 


57:17.2 


+49:17:08 


17.0 [T88] 


262 


12 


38 


-18.77 


7 





776.7 + 0.7 


121+ 20 


39.5 ± 0.5 


25 + 4 


140 


NGC 4519 


12 


33:30.5 


+08:39:16 


19.6 [T08] 


190 


12 


15 


-19.31 


7 





1218.1 ± 1.0 


112+7 


355 ± 2 


42 + 3 


218 


NGC 4561 


12 


36:08.6 


+ 19:19:26 


12.3 [T88] 


90.8 


12 


82 


-17.63 


8 





1402.2 ± 0.9 


57+4 


227 ± 8 


34 + 4 


171 


NGC 4713 


12 


49:58.1 


+05:18:39 


14.9 [T08] 


162 


11 


85 


-19.02 


7 





654.5 ± 0.5 


110. 9+1 s 


274.0 ± 0.6 


45.2 ± 1.2 


176 


NGC 4942 


13 


04:19.2 


-07:39:00 


28.5 [T88] 


112 


13 


27 


-19.00 


7 





1741 ± 2 


124t 6 5 ' 


137.3 ± 0.8 


37 + 2 


177 


NGC 5964 


15 


37:36.3 


+05:58:28 


24.7 [T88] 


250 


12 


28 


-19.68 


7 





1447.1 ± 1.2 


168+;* 


136.7+ 1.2 


32 + 3 


208 


NGC 6509 


17 


59:24.9 


+06:17:12 


28.2 [T88] 


95.1 


12 


12 


-20.13 


7 





1811.0 ± 0.4 


153+ 12 


280.8 ± 1.1 


41 + 4 


266 


IC 1291 


18 


33:51.5 


+49:16:45 


31.5 [T88] 


109 


13 


28 


-19.21 


8 





1951.0 ± 1.1 


189111 


131 ± 2 


28 + 3 


209 



Note. — Column 1: Object name; Column 2 and 3: Right ascension and declination (J2000.0) fr om Ide Va^o^uleurB et al.l iH)91|) ; Colu mn 4: Distance and 
distance reference. Distances are derived using th e Tully-Fisher relation, except for NGC 3906. T08: iTullv et alj teOOSft. T 88: iTullvl dl988D . and the NGC 3906 
distance is from the |de Vaucouleurs et alj dl991fl heliocentric velocity, correc ted for Virgo infall usin^TMould^t**alT ~ ?2000|) and using H Q = 71 kms -1 Mpc -1 . 
Column 5: Major isophotal diameter at 25 mag arcsec - 2 in the B band, fromlde V auco uleurs et al.l j*199llv Column 6: Apparent blue magnitude, corrected for 
Galactic and internal extinction and rcdshift. Values arc from lde Vaucouleurs et alj U991D . except for NGC 4942 and PGC 6667, which are from [Doyle et alj 
(2005) and arc only corrected for Galactic extinction. Column 7: Absolute b lue m agnitude, calculated from the apparent magnitude in Column 6 and the distance 
in Column 4. Column 8: Morphological type from dc Vaucouleurs ct al. (1991). Column 9: Systemic velocity, corrected to the heliocentric reference frame. 
Column 10: Circular velocity. Columns 9 and 10 were derived from the VLA H I rotation curve analysis in Paper I. Column 11: Position angle of the major axis 
(degrees N to E to receding side). Column 12: Inclination. Columns 11 and 12 were derived from a combination of the rotation curve analyses on the H I data 
and ellipse fits of the IRAC 3.6 f-Lm. data. Columns 9-12 were originally presented in Table 6 of Paper I. Column 13: Width of the HI line at 20% of the peak flux 
density, corrected for the spectral resolution but not turbulent broadening. The uncertainty in W20 is 5 kms" 1 for all objects except UGC 6446 and NGC 3906, 
where the uncertainty is Skins" 1 and 10 kms -1 , respectively. Column 13 was originally presented in Table 4 of Paper I. 



to 
to 

TABLE 2 

Measured and Literature Properties 



Source 
(1) 


RA 
(hh:mm:ss.s) 
(2) 


DEC 
(dd:mm:ss) 
(3) 


(10~ 14 ergs -1 cm -2 ) 
(4) 


fpAH 

(mjy) 
(5) 


Chi) 

( Jy beam 1 km s 1 ) 
(6) 


Bmaj 

(arcscc) 
(7) 


-Bmin 

(arcscc) 
(8) 


^co 
(K km s~ 

(9) 


) 


F 4.5 

(mJy) 
(10) 


F 4.5 

(mJy) 
(11) 




J 

(mag) 
(12) 


K s 
(mag) 
(13) 


NGC 0337 


00:59:49.9 


-07:34:44 


10 9 + 005 


73 + 8 


1.28 ± 


17 


25.77 


15 


.06 


4.7 ± 0. 


3 


11.7 ± 1.2 


64 + 6 


9.876 ± 0.025 


9.059 ± 0.045 


PGC 3853 


01:05:04.9 


-06: 12:45 


(I QS + fl Q4 


6.9 ± 0.8 


0.36 ± 0. 


05 


25.33 


15 


.47 


1.39 ± 0. 


18 


3.4 ± 0.3 


44 + 4 


10 


.031 ± 0.037 


9.280 ± 0.080 


pnr 1 6667 

JT VJU UUU 1 


01 -4Q- 1 n 3 


1 0-03-40 


o 97 _|_ n no 


8 5 + 09 


0.90 ± 0. 


12 


26 00 


21 


.00 


71 + 


13 


2 4 + 02 


15.2 + 15 


11 


.678 ± 0.037 


1 n Q51 -i- n 080 




02-1 4- ^6 & 


-20-12-44 


n sc _|_ n no 


5 3 + 06 


0.59 ± 0. 


.08 


27 78 


12 


.01 


< 4 




2 5 + 03 


12 8 + 1 3 








TIGC 1862 


02:24:24.8 


-02:09:44 


-i Q7 J. n no 
Iil/ i _i_ yj . u z- 


4.2 ± 0.5 


0.24 ± 


03 


23.97 


15 


.32 


0.60 ± 0. 


10 


2.2 ± 0.2 


12.4 ± 1.2 


11 


.927 + 0.023 


11 1 77 + 054 
11.11 1 _1_ y.f .Kj'.j'rt 


rso 4i 8-0008 


03-31-30 7 


-30:12:48 


1 04 _L n no 


6.7 ± 0.8 


0.92 ± 0. 


12 


32.00 


21 


.00 


< 0.9 




2.7 ± 0.3 


6.7 ± 0.7 


12 


.752 ± 0.049 


12.169 ± 0.126 


ESO 555-0027 


06-03-36 8 


-90-3Q-1 5 


1 49 + D 05 

J. . '-t _1_ k J . U 


9.3 ± 1.0 


0.44 ± 0. 


06 


29.00 


21 


.00 






2.4 ± 0.2 


16.1 ± 1.6 


11 


.984 ± 0.045 


1 1 271 + 1 03 


NOf! 2805 


09-20-20 3 


+ 64-06-1 1 


w _|_ n m 7 


15.5 ± 1.7 


0.34 ± 


04 


21.00 


21 


.00 


2.34 ± 0. 


10 


4.2 ± 0.4 


41 + 4 


10 


.827 ± 0.026 


in 1 1 7 _|_ n nAP. 
iu*ii 1 _i_ y.f . yj t: yj 


ESO 501-0023 


1 0-35-23 3 


-24:45:15 


f) 76 + 03 

O. 1 U _1_ U . KJ O 


1.04 ± 0.14 


0.134 ± 


017 


30.30 


12 


.51 


< 0.6 




90 + 09 


9.9 ± 1.0 








UGC 6446 


11:26:40.4 


+53:44:48 


1.19 ± 0.02 


0.44 ± 0.09 


0.43 ± 


06 


21.00 


21 


.00 


< 0.4 




0.97 + 0.10 


9.5 ± 1.0 








NGC 3794 


11:40:54.3 


+56:12:07 


2.37 ± 0.016 


7.3 ± 0.9 


0.49 ± 0. 


06 


21.00 


21 


.00 


0.72 ± 0. 


15 


3.5 ± 0.4 


13.0 ± 1.3 


11 


.871 ± 0.039 


11.005 ± 0.072 


NGC 3906 


11:49:39.9 


+48:25:32 


2.25 ± 0.02 


7.7 + 0.9 


0.25 ± 


03 


21.00 


21 


.00 


0.37 ± 0. 


12 


2.8 ± 0.3 


14.0 ± 1.4 


11 


.788 ± 0.045 


11.007 + 0.067 


UGC 6930 


11:57:17.4 


+49:16:58 


0.76 ± 0.03 


6.8 ± 0.8 


0.25 ± 0. 


.03 


21.00 


21 


.00 


1.31 ± 0. 


19 


2.6 ± 0.3 


29 + 3 


11 


.752 ± 0.037 


11.153 ± 0.068 


NGC 4519 


12:33:30.3 


+08:39:18 


3.0 ± 0.03 


44 + 5 


0.99 ± 0. 


13 


51.91 


18 


.77 


2.9 ± 0. 


3 


6.9 ± 0.7 


36 + 4 


10 


.499 ± 0.037 


9.548 ± 0.059 


NGC 4561 


12:36:08.2 


+ 19:19:22 


4.55 ± 0.03 


7.9 ± 1.0 


1.02 ± 


13 


21.00 


21 


.00 


0.68 ± 0. 


12 


4.1 ± 0.4 


14.8 ± 1.5 


11 


.480 ± 0.046 


10.617 + 0.074 


NGC 4713 


12:49:58.0 


+05:18:41 


5.35 ± 0.03 


35 + 4 


0.69 ± 0. 


09 


21.00 


21 


.00 


4.5 ± 0. 


2 


7.1 ± 0.7 


44 + 4 


10 


.367 + 0.026 


9.737 ± 0.053 


NGC 4942 


13:04:19.1 


-07:38:58 


3.55 ± 0.04 


10.1 ± 1.1 


0.49 ± 0. 


.06 


22.00 


21 


.00 


1.4 ± 0. 


2 


3.4 ± 0.3 


12.9 ± 1.3 


11 


.617 + 0.058 


10.428 ± 0.092 


NGC 5964 


15:37:36.2 


+05:58:27 


2.11 ± 0.05 


9.0 ± 1.0 


0.30 ± 0. 


.04 


21.00 


21 


.00 


0.89 ± 0. 


17 


3.4 ± 0.3 


40 + 4 


12 


.384 ± 0.058 


11.794 ± 0.113 


NGC 6509 


17:59:25.2 


+06:17:11 


4.62 ± 0.03 


41 + 5 


0.50 ± 0. 


.07 


21.00 


21 


.00 


6.0 ± 0. 


2 


10.6 ± 1.1 


42 + 4 


10 


.284 + 0.024 


9.663 ± 0.043 


IC 1291 


18:33:52.5 


+49:16:42 


6.02 ± 0.02 


9.3 ± 1.0 


0.70 ± 0. 


09 


21.00 


21 


.00 


0.37 ± 0. 


08 


2.4 ± 0.2 


10.2 ± 1.0 


13 


.125 ± 0.052 


12.689 ± 0.152 



Note. — Column 1: Object name. Column 2 and 3: RA and DEC (J2000.0) of the pointing center of the IRAM 30m CO(1-0) observations. All measurements within a 2l"-diameter circular aperture 
are centered on these coordinates. Column 4: Ha flux within a 2l"-diameter circular aperture, corrected for Galactic extinction and N II emission. We quote only measurement uncertainties here. Column 
5: PAH flux density within a 2l"-diameter circular aperture. Column 6: Average integrated H I line intensity within a 2l"-diameter circular aperture, measured from image convolved to have beam major 
and minor axes given in Column 7 and 8. Column 7 and 8: Beam major and minor axes of image used for HI line intensity measurement. Column 9: Integrated CO(1-0) line intensity. Non-detections arc 
quoted as 3c upper limits. Column 10: 4.5/j.m flux density within a 21 -diameter circular aperture. Column 11: Total 4.5/im flux density. Column 12: Total J-band magnitude from 2MASS, corrected 
for Galactic extinction. Column 13: Total i\g-band magnitude from 2MASS, corrected for Galactic extinction. 
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TABLE 3 
Ha Observations 



Galaxy 


*663ni)15 (min) 


Run(s) 


*693nbl5 (min) 


Run(s) 


(i) 


(2) 


(3) 


(4) 


(5) 


NGC0337 


2400 


3 


2400 


3 


PGC3853 


1800 


4 


1800 


1 


PGC6667 


4500 


1,3 


4500 


1,3 


ESO544-G030 


5400 


1,3 


5400 


1,3 


UGC1862 


4500 


1,3 


5400 


1,3 


ESO418-G008 


2700 


1 


2700 


1 


ESO555-G027 


2700 


1 


2700 


1 


NGC2805 


2400 


3 


2400 


3 


ESO501-G023 


3600 


1 


3600 


1 


UGC6446 


3900 


2 


4200 


2 


NGC3794 


6000 


1,3 


6000 


1,3 


NGC3906 


3200 


1,2,3 


3000 


1,2 


UGC6930 


3600 


2 


3600 


2 


NGC4519 


3600 


2 


3600 


2 


NGC4561 


3600 


1 


3600 


1 


NGC4713 


2400 


1,2 


2400 


1.2 


NGC4942 


4500 


1.2 


3600 


2 


NGC5964 


3600 


2 


3600 


2 


NGC6509 


4800 


2 


3600 


2 


IC1291 


3600 


2 


3600 


2 



Note. — Column 1: Object name. Column 2: Total exposure time 
in the 663nbl5 filter. Column 3: Observing run code. Column 4: 
Total exposure time in the 693nbl5 filter. Column 5: Observing run 
code. Observing run codes are: (1) Jan 2007; (2) May/June 2007; (3) 
Nov 2007; (4) Jan 2008. 



TABLE 4 
Derived Properties 



Source 

(i) 


v 

(2) 


log S H i 

(3) 


log S H 2 

(4) 


logS SFR 
(5) 


logS, 
(6) 


log M, 
(7) 


12 + log(0/ff) 
(8) 


U 
(9) 


K 

(10) 


Qgas 
(11) 


Q, 
(12) 


Qgas+stara 

(13) 


log ft 
(14) 


NGC 0337 


2 


.11 


1 


.33 


1.18 


-1.67 


2 


.49 


9.9 


8.93 


1.7 


0.17 


2.7 


2.9 


2.0 


5.32 


PGC 3853 


1 


.16 





.78 


0.67 


-2.68 


1 


.88 


9.1 


8.68 


1.6 


0.22 


12 


8 


5 


4.42 


PGC 6667 


2 


.51 


1 


.09 


0.42 


-2.44 


1 


.73 


9.3 


8.75 


2.4 


0.12 


4.8 


6 


3.2 


4.50 


ESO 544-G030 


1 


.42 


1 


.02 


< 0.09 


-2.83 


1 


.51 


8.6 


8.46 


1.1 


0.12 


0.0 


5 


2.9 


4.39 


UGC 1862 


2 


.27 





.62 


0.29 


-2.70 


1 


.67 


9.2 


8.70 


1.5 


0.07 


7 


3.0 


2.3 


4.04 


ESO 418-G008 


2 


.40 





.84 


< 0.34 


-2.72 


1 


.40 


8.7 


8.50 


1.1 


0.07 


4.5 


3 


2.0 


4.21 


ESO 555-G027 


2 


.47 





.78 




-2.44 


1 


.77 


9.3 


8.74 


1.9 


0.12 




5 






NGC 2805 


2 


.85 





.74 


0.92 


-2.28 


1 


.93 


9.8 


8.91 


2.3 


0.041 


1.8 


1.6 


1.1 


4.52 


ESO 501-G023 





.71 





.40 


< 0.32 


-3.16 


1 


.01 


7.8 


8.00 


0.72 


0.12 


16 


8 


5 


3.77 


UGC 6446 


1 


.83 





.73 


< -0.008 


-3.21 


1 


.01 


8.6 


8.47 


1.1 


0.07 


6.5 


5 


2.9 


3.88 


NGC 3794 


1 


.96 





.75 


0.27 


-2.64 


1 


.93 


0.2 


8.72 


1.5 


0.15 


12 


5 


3.9 


4.26 


NGC 3906 


1 


.86 





.69 


0.20 


-2.42 


1 


.94 


9.1 


8.67 


0.93 


0.09 


8 


2.1 


1.8 


4.27 


UGC 6930 


1 


.73 





.66 


0.72 


-2.63 


1 


.62 


9.1 


8.66 


2.1 


0.18 


11 


10 


6 


4.24 


NGC 4519 


2 


.00 





.82 


0.99 


-1.93 


2 


.46 


9.8 


8.90 


2.0 


0.22 


8 


4 


3.7 


4.83 


NGC 4561 


1 


.25 


1 


.23 


0.40 


-2.32 


2 


.15 


8.9 


8.58 


0.69 


0.14 


4.2 


2.3 


1.7 


4.98 


NGC 4713 


1 


.52 





.99 


1.15 


-1.98 


1 


.99 


9.2 


8.70 


1.4 


0.25 


6.3 


7 


4.0 


4.94 


NGC 4942 


2 


.90 





.88 


0.69 


-2.35 


2 


.52 


10.0 


8.96 


1.8 


0.12 


5.8 


2.1 


1.9 


4.72 


NGC 5964 


2 


.52 





.71 


0.53 


-2.44 


1 


.70 


9.5 


8.82 


3.4 


0.055 


3.9 


3 


2.2 


4.11 


NGC 6509 


2 


.87 





.89 


1.31 


-1.92 


2 


.18 


9.7 


8.88 


1.7 


0.22 


4.7 


6 


3.3 


5.06 


IC 1291 


3 


.21 


1 


.10 


0.17 


-2.19 


1 


.35 


8.9 


8.60 


2.1 


0.18 


8 


13 


5 


4.37 



Note. — Column 1: Object name. Column 2: Physical size of 21" (kpc). All surface density measurements arc within a 21 // -diamctcr aperture. 
Column 3: HI mass surface density (Mq pc" 2 ). The typical uncertainty is 0.06 dex. Column 4: H2 mass surface density (Mq pc~ 2 ), computed with 
Xqo — 2. 8 X 10 20 cm" 2 ( K km s" 1 )" 1 . The typical uncertainty is 0.07dex. Column 5: SFR surface density, computed assuming the Kroupa-typc 
IMF from lCalzetti et al.| J2007I) (M yr _1 kpc" 2 ). Wc assume a typical uncertainty of 0.2 dex. Column 6: Stellar mass surface density {Mq pc 2 ). 
We assign a a typical uncertainty of 0.2 dex. Column 7: Total stellar mass (M q). We assign typical uncertainty of 0.3 dex. Column 8: Oxygen 
abundance derived from the total stellar mass and the mass-metallicity relation of Trcmonti ct al. ( 2004). Wc derive a typical statistical uncertainty 
of 0.11 dex, but note that strong-line metallicity methods return 12 -f \og{0 / H) values as disparate as 0.5 dex. Column 9: Stellar scale length 
(kpc). Wc estimate the uncertainty to be 20%. Column 10: Epicyclic frequency (km/s/pc), evaluated at 5.25 . We estimate the uncertainty to 
be 30%. Column 11: Gas stability parameter. The typical uncertainty is 30%. Column 12: Stellar stability parameter. The typical uncertainty is 
60%. Column 13: Combined gas and stellar stability parameter. The typical uncertainty is 40%. Column 14: Mid-plane pressure (Kern -3 ). The 
typical uncertainty is 0.17 dex. 



